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6. FOLGE * BAND 15, HEFT 1 1954 


Unter welchen Bedingungen ist der Amplitudenverlauf 
einer seitlich begrenzten Welle komplex? 


Von Kurt Artmann 


Inhaltsübersicht 


Die spezielle Form des von 0 auf 1 ansteigenden Amplitudenverlaufs 
A(l,b) einer seitlich begrenzten (Licht-)Welle A(l,b)ei*! kann in allen 
praktischen Fällen nicht willkürlich (reell) vorgegeben werden, sondern muß 
in Gestalt des (wesentlich komplexen) Fresnelschen Integrales dargestellt 
werden. 


Bezeichne / bzw. b die Koordinate in bzw. senkrecht zur Fortschreitungs- 
richtung der Welle, so wird in vielen Untersuchungen die Wellengleichung 
+ + k? y = 0 (k = 2 n/A) durch den Ansatz 


y= A (I, b)eikl (1) 


gelöst. In einer früheren Arbeit!) zeigte der Verfasser, daß die Wellenglei- 
chung für alle J, b nur dann durch eine reelle Amplitudenfunktion A (l, b) 
befriedigt werden kann, wenn A = const. b + const. Aus diesem Grunde 
hatte der Verfasser damals Einwände gegen eine Arbeit von Cl. Schaefer, 
Pich?) sowie v. Fragstein, Schaefer?) erhoben, in denen mit einem 
reellen A (l, b) operiert wurde. 

Diese Einwände sehen v. Fragstein und Schaefer‘) als ungerecht- 
fertigt an und erheben außerdem Vorwürfe gegen die Richtigkeit von zwei 
Arbeiten des Verfassers. Diese Gegenargumente — soweit sie nicht Prioritäts- 
ansprüche darstellen — fußen auf folgenden beiden mathematischen Be- 
hauptungen, von denen die eine (A) sicher unrichtig ist. Die zweite Behauptung 
(B) ist zwar wörtlich genommen richtig. Aber ihre Formulierung ist so gewählt, 
daß sie das Gegenteil von dem ausdrückt, was wirklich stets in der Praxis 
vorliegt (nämlich komplexes 4). 


A. Die Superposition zweier ebener Wellen mit gleicher Wellenlänge 


4 =2a/k aber etwas verschiedenen Fortschreitungsrichtungen läßt sich in 
Strenge darstellen in der Form 


y=eikl.F(b) (2) 


1) K. Artmann, Ann. Physik (6) 8, 285 (1961); im folgenden als I zitiert. 

2) Cl. Schaefer u. R. Pich, Ann. Physik (5) 30, 245 (1937); im folgenden als II 
zitiert. 
- 3) Cl. Schaefer u. C. v. Fragstein, Ann. Physik (6) 6, 39 (1949); im folgenden als 

zitiert. 

4) C. v..Fragstein u. Cl. Schaefer, Ann. Physik (6) 12, 84 (1953), im folgenden 
als IV zitiert. 
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mit reellem F (6). (IV, S. 87 oben). Hierin bedeutet J die Koordinate in 
Richtung der Winkelhalbierenden der beiden ursprünglichen Wellen; 5 ist 
die Koordinate senkrecht hierzu. 
B. Unter der Voraussetzung 
RA | dA 
| <|A| (3) 
läßt sich die Wellengleichung näherungsweise durch reelles A (1, b) erfüllen 


(IV, S. 85 unten). 
Zu A. Die Superposition zweier ebener Wellen mit den Fortschreitungs- 


richtungen y + 4/2: 
y= eiklzsin 4/2) -zc08 (g-4/2)} + 
ergibt nach Einführung der Abkürzungen 
l= axsing—zcosy, b= xcosg + zsing 
nach elementarer Zwischenrechnung ohne irgendwelche Streichungen 
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4 
y= F(l,b) = 2-cos(kbsin 5) (4) 


Wegen der Exponentialfunktion ist F nicht reell, also die Behauptung A 
widerlegt. 3 

Entwickelt man bei kleinem A den cos 4/2 der Exponentialfunktion von 
(4), so nimmt diese die Form e~‘*!4*/8 an. In IV wird die Exponential- 
funktion wegen A< 1 von v. Fragstein und Schaefer durch den Wert 1 
ersetzt, wodurch dann F reell ausfällt. Hierbei ist nicht beachtet, daß diese 
Vereinfachung nur innerhalb eines Gebietes der Lineardimension | < 8 2/(k A?) 
erlaubt ist. Da der Verfasser im § 2 von V°) zwei ebene Wellen super- 
ponierte, welche den gesamten Raum durchziehen (-oo <1< + oo), 
so ändert sich die Phase von F (l, b) wesentlich im gesamten Raume. Meinen 
damaligen Betrachtungen lag also ein komplexer Verlauf F (1, b) [=A (I, b)] 
zugrunde. Damit entfallen die in IV, S. 87, 1. Abschnitt, gemachten Vorwürfe 
gegen die Richtigkeit meiner Betrachtungen in V. 

Da die Phase von F nur von der Längskoordinate / abhängt, und außerdem 
innerhalb hinreichend kleiner Raumgebiete praktisch konstant ist, so braucht 
diese Phasenveränderlichkeit bei den meisten Betrachtungen nicht mit- 
berücksichtigt zu werden. Da aber die v. Fragstein-Schaeferschen Ein- 
wände gegen die Richtigkeit meiner früheren Überlegungen allein auf dem 
(falschen) Argument aufgebaut sind, daß F in aller Strenge reell sei, so 


‚muß bei der vorliegenden Entgegnung ausdrücklich auf diese Phasenver- 


änderlichkeit hingewiesen werden. 

Hinzu kommt, daß die auf der Superposition von nur zwei ebenen Wellen 
basierende elementare Theorie im $2 von V nur dazu dienen sollte, auf denk- 
bar kürzestem Wege den Wert für die Strahlversetzung beim Versuch von 
Goos und Hänchen®) verständlich zu machen. Meine strenge Theorie 
(Teil II von V), bei welcher unendlich viele ebene Wellen superponiert 
wurden, operierte mit dem Fresnelschen Integral als Amplitudenfunktion 


5) K. Artmann, Ann. Physik (6) 2, 87 (1948); im folgenden als V zitiert. 
°) F. Goos, H. Hänchen, Ann. Physik (6) 1, 333 (1947). 
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A (l,b), die wesentlich komplex ist (vgl. die folgende Anm. 7). Leider er- 
wähnten v. Fragstein und Schaefer in IV, 8. 87 oben, diesen Umstand 
nicht. 

Zu B. Dem Ansatz (1) entspricht es, sich auf Amplitudenverläufe A (J, b) 
zu beschränken, die langsam auf der Strecke einer Wellenlänge veränderlich 
sind, d.h. 

aA aA aA 0A 

9 

Den schnell oszillierenden Teil von y denkt man sich mit in den Faktor e‘*! 

in (1) aufgenommen. Da y der Wellengleichung genügt, so befriedigt A die 
Differentialgleichung 


aA 
+ = 0. (6) 


Wegen (5) kann der erste Term in (6) gegen den dritten gestrichen werden, 
so daB 
aA 
+ 2ik 0. (7) 
Bezeichnet A (0, b) den vorgegebenen Wert in einer Anfangsebene / = 0, 
so läßt sich die Lösung von (7) bekanntlich darstellen in der Form 


A (I, b) = eine. Vers A(0,b)e 2 db. (8) 


Im Spezialfall, daß A(0,b’) der Dirichletfaktor A (0,b') =1 für b'> 0, 
A (0, b’) = 0 für b’ < 0, ist (8) ersichtlich identisch mit dem Fresnelschen 
Integral 


. _ byk 
Praktisch ist A (l,b) immer dann gleich dem Fresnelschen Integral (9), 
wenn A (0, b') rasch ansteigt innerhalb des Gebietes (b — b’)? k/(2 1) S 2/8, 
innerhalb dessen der 2. Faktor in (8) sich nur wenig ändert. Bezeichnet b, 
den Bereich des wesentlichen Anstiegs der Amplitudenfunktion A (0, 5) in 
der Anfangsebene != 0, so ist die Amplitudenfunktion A(l,b) im 
Abstande 


Adl,b) = 


LRA (10) 


von der Anfangsebene praktisch immer gleich dem Fresnelschen Integral 
unabhängig von der speziellen Form des Amplitudenanstiegs 
A(0,b) auf der Anfangsebene 1=0. Da das Fresnelsche Integral 
wesentlich komplex’) ist, so ist auch A (l,b) wesentlich komplex, 
sofern (10) erfüllt ist. 


?) Man erkennt dies am einfachsten, wenn man die im Bereich |&] S 0,7 (d. h. in Nähe 
der geometrischen optischen Schattengrenze) gültige Näherungsdarstellung von (9) 


benutzt: A = + en eg, 
2 Ya 
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Da die Amplitude A (l,b) hiernach nach Durchlaufen der letzten Blende 
immer bestrebt ist, die Form eines Fresnelschen Integrales anzunehmen, so 
ist bei einer Versuchsanordnung vom Typ Goos-Hänchen, bei welcher der 
Strahl nach Durchlaufen der letzten Blende wegen der vielen Zickzack- 
reflexionen einen langen optischen Weg / zurücklegt, die Bedingung (10) 
immer erfüllt, also A (l,b) gleich dem komplexen Fresnelschen Integral. 
Mit dem Werte (9) geht die v. Fragstein-Schaefersche Voraussetzung 
(3) im allein interessierenden Bereich |£| s 0,7 über in die Forderung kl>1, 
die sicherlich erfüllt ist. Damit ist die Behauptung B widerlegt; und somit 
entfallen die in IV, S. 86, 2. Absatz, gegen mich vorgebrachten Vorwürfe. 

Um die v. Fragstein-Schaefersche Voraussetzung zu erfüllen, daß 
A (l,b) auf der Plattenebene reell ist, müßte durch passende Versuchsan- 
ordnung erreicht sein, daß nach Durchlaufen der letzten Blende die Ampli- 
tudenfunktion A (0,5) in der Anfangsebene / = 0 (praktisch) rein reell ist 
mit einem Anstiegsbereich ux Vik (l ist der optische Weg zwischen der 
Anfangsebene und der Aufnahmeplatte). Beim Versuch von Goos-Hänchen 
war 2» 10cm. Mit k = 10° em”! müßte an 10? em sein. Bedenkt man, 
daß im Abstande einiger weniger Wellenlängen hinter jeder metallischen 
einbackigen Blende schon ein (komplexer) Fresnel-Verlauf A (l, b) vorliegt, 
so dürfte es nur unter völlig anomalen Versuchsbedingungen möglich sein, 
auf der Aufnahmeplatte einen reellen Amplitudenverlauf A (l,b) herzustellen, 
der von 0 auf 1 ansteigt. Außerdem wirkt der soeben errechnete große Wert 
von b, der Forderung entgegen, daß die Flankensteilheit von A möglichst 
klein sein muß, um den eigentlichen Effekt der Strahlversetzung möglichst 
genau auszumessen. 

Die Diskrepanzen, die in IV zwischen meinen Ergebnissen und denen von 
Wolter®) aufgezeigt wurden, beruhen auf dem Umstande, daß ich gemäß 
dem in II und III aufgestellten Programm einen Amplitudenverlauf A (I, b) 
behandelte, der von 0 auf 1 ansteigt, während bei der Wolterschen Versuchs- 
anordnung der unter A diskutierte Amplitudenverlauf (4) zugrunde lag, der 
den gesamten Raum durchzieht. Da A (l,b) im ersten Fall wesentlich, 
im zweiten Fall nur unwesentlich komplex ist, so ist es offenkundig, daß 
man bei einer Diskussion der Realitätsverhältnisse von A (l, b) die Wolter- 
sche Versuchsanordnung nicht ohne weiteres mit der von Goos-Hänchen 

vergleichen darf, wie es in IV geschieht. 

Es bleibt somit das Ergebnis bestehen, daß von den bisher praktisch 
durchgerechneten Fällen in der Schaefer-Pichschen Theorie nur bei der 
Wolterschen Versuchsanordnung mit reellem A (l, b) operiert werden darf. 
Diese Tatsache schränkt den Wert der Schaefer-Pich-v.-Fragsteinschen 
Theorie insofern erheblich ein, als sich der Woltersche Fall, bei dem gemäß A 
nur zwei ebene Wellen superponiert werden, in fast trivial kurzer Weise 
durchrechnen läßt, was bei Schaefer-Pich sicher nicht der Fallist. Der Sinn 
der Rechnungen von II und II kann doch nur darin bestehen, auf Probleme 
angewendet zu werden, bei denen nicht zwei, sondern unendlich viele 
ebene Wellen superponiert werden müssen, wie es bei dem in II und III auf- 
gezeichneten Amplitudenverlauf A der Fall ist. 


8) H. Wolter, Z. Naturforschg. 5a, 143 (1950). 
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Nach dem Vorangegangenem ist es dann aber nicht — wie v. Fragstein 
und Schaefer meinen — möglich, allein durch vereinfachte Annahmen über 
den Amplitudenverlauf A (I, b) die ,, Beugung zu eliminieren“ (IV, S. 87 unten), 
wenn die Rechnungen auf Experimente mit einer seitlich begrenzten Welle 
angewendet werden sollen. Auf die Schrödingersche Wellenmechanik 
übertragen, würde diese „Elimination der Beugung‘‘ darauf hinauskommen, 
das Zerlaufen eines Wellenpaketes im Laufe der Zeit völlig zu ignorieren. Denn 
die zeitabhängige Schrödinger-Gleichung, eines in einer (b-)Richtung 
beweglichen Wellenpakets, ist mit (7) identisch, sofern in (7) die Längskoor- 
dinate | durch A kt/m ersetzt wird (m = Teilchenmasse, # = Plancksche 
Konstante). Bei einer üblichen (optischen) Versuchsanordnung, bei welcher 
A (l, b) praktisch gleich dem Fresnelschen Integral ist, bedeutet eine solche 
von Schaefer-Pich-v.-Fragstein vorgenommene „Elimination der Beu- 
gung durch Annahme eines reellen Amplitudenverlaufs A (I, b)‘‘ nichts weiter 
als einen wesentlichen Verstoß gegen die Wellengleichung mit den 
von mir in I aufgezeigten Konsequenzen, die ich nach wie vor aufrechterhalte. 


Hamburg, Institut für Theoretische Physik der Universität. 


Bei der Redaktion eingegangen am 1. März 1954. 
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Effect of Acoustic Impedance and Viscosity of Gases 
on the Electrical Constants of Quartz 


By S. Parthasarathy and V. Narasimhan 
With 6 Figures 


Abstract 


The equivalent resistance, capacitance and the Q factor of a quartz crystal 
vibrating at its fundamental frequency in several gaseous media were measured 
and the results compared with the acoustic impedance (P C) of the media. 
While the resistance varies linearly with P C the reciprocal capacitance does 
not. The Q factor seems to vary linearly with viscosity within the small range 
of viscosities. 


1. Introduction 


The electrical equivalent circuit of the piezo-electric transducer has been 
worked out by Cady"), Van Dyke?) and Dye’). It has been shown that the 
quartz element can be assumed to be made up of a resistance, an inductance 
and a capacity all in series. Hubbard‘) generalised the theory of the equi- 
valent network and deduced the conditions under which one could estimate 
the reflection and attenuation coefficients when the quartz element was vibra- 
ting in a medium. Fox 5) has applied Hubbard’s theory for the calculation 
of absorption coefficients of ultrasound in water. Cady!) has shown that 
when the quartz crystal vibrates in one of its 
resonant frequencies, the current passing 
through it drops considerably and under 
ideal conditions may reach even near zero 
values. 

Hubbard considers the effect of the 
coupled fluid column on the electrical equi- 
valents of the piezo crystal. He has shown 
that the coupled fluid column modifies the 
equivalent resistance and the piezoelectric capacity of the equivalent circuit 
of the quartz resonator. Referring to fig. 1, Hubbard obtains the expression 


o2 


2 2 
1) W. G. Cady, Proc. Inst. Radio Engineers 10, 83 (1922). 
2) K. S. van Dyke, Proc. Inst. Radio Engineers 16, 742 (1928). 
») D. W. Dye, Proc. physic. Soc. London 88, 39 (1926). 
4) J.C. Hubbard, Physic. Rev. 88, 1011 (1931); 41, 523 (1932); 46, 525 (1934). 
5) F. E. Fox, Physic, Rev. 52, 973 (1937). 
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where L,, C,, R,, are the inductance, capacity and resistance of the pick up 
electrical resonant circuit K,, K’, R’ and L are the corresponding equivalent 
electrical constants of the quartz resonator with the coupled fluid column K, 
being the dielectric capacity of the resonator and K’ the modified piezo- 
electric capacity and 

P =1-(4 +h) 

=I-LK’' 

®, = R,(C, + Ky) 

®, = R' K'w 

o =i/I, 
and i is the current in the Z, R, branch. Jy = 4° is the maximum current in 

1 

the resonant circuit when the crystal is disconnected E, being the maximum 
voltage induced in L,, R’ and K’ include the modifications produced by the 
fluid column and 

R’=R+ABPcS (2) 


(3) 


where R is the equivalent resistance of the piezo crystal near its response 
frequency, K its piezo-electric capacity, A the effective area of the crystal 
exposed to the fluid column, B a piezo-electric coefficient, P the density of 
the fluid, C the velocity of sound in it and @ is 2 // times the frequency of the 
voltage impressed in the coil. 

S and T are further defined so as to include terms involving the attenua- 
tion factor of the particle velocity in the fluid, the coefficient of reflection at 
the distant boundary and the length of the fluid column (Hubbard). When 
the sound wave is absorbed by the distant boundary the terms S and T 
become each equal to unity and equations (2) and (3) reduce to the simple 
form 

R=R+ABPC (4) 


K' 


=ZtABPC. (6) 


Equations (4) and (5) are of special interest, since the equivalent resistance 
and reciprocal capacitance of a piezo crystal vibrating in different media 
should vary directly with the acoustic impedance P C of the medium and so 
the following experiments were undertaken to verify how far these relations 
could be assumed to represent conditions actually obtaining when a piezo 
crystal vibrates in several gaseous media of varying P C. 


2. Experimental 


Several preliminary investigations were conducted to obtain a good cre- 
vasse but the one which was eventually successful will be described. 


Ss 
stal 
ıred 
dia. 
loes 
Inge 

1 1 i 
een 
the | 
nee 
qui - 
late 
ora- 
ion | 
hat 
its 
ing 
der 
ero 
> 
ui- 
wn 
the 
uit | 
ion | 
(1) 


Annalen der Physik. 6. Folge. Band 15. 1954 


A Marconi TF 144 G type signal generator excited the crystal in its 
fundamental frequency. This signal generator has been claimed by its makers 
to be of substandard accuracy. The signal from this generator was fed into 
a General Radio unit amplifier. The output of this amplifier was fed to the 
coupling coil L, such that the amplifier delivered power into the loosely 
coupled L, C, circuit tuned to the crystal frequency, fig. 2. This method 
is the one generally adopted for all ultrasonic interferrometry. L, is a coil of 
low inductance and was a standard 25 wH coil supplied by Marconi’s. The 
vacuum tube voltmeter capable of measuring voltages as low as 15 millivolts 
was used to measure the voltages across the crystal. The readings of this 
instrument were compared with those of the standard voltmeters available 
in this laboratory. For 


the measurement of ex- 

SIGNAL med tremely low voltages that 

GEN. ‚ame | ı3 c, oceur at the bottom of a 

crevasse, asensitiveshunt- 

ed galvanometer was con- 

X- CRYSTAL nected across the micro- 

Fig. 2 ammeter of the vacuum 

tube voltmeter through a 

changeover commutator and readings could be taken either on the voltmeter 
or galvanometer at will. 


It was the aim of the experiment to oscillate the quartz crystal in several 
gaseous media and arrangements were made to evacuate and fill the gas 
container. The set up is described in fig. (3). These gases were dried and 
freed from carbon dioxide before letting into the chamber. The pressure was 
read on a barometer incorporated in the system and the temperature was 

r measured by a thermometer with 
> 79 pump, ts bulb close to the quartz crystal. 
DRYING TUBES Of the several X-cut quartz 
crystals available here, one showed 
an excellent crevasse. It is well- 
that X-cut crystals always 
D.Dewarrıssa one or two secondary dips 
M.MANOMETER associated with them and the one 
8-GAROMETER which has a dip on one side of the 

Fig. 3 principal dip had to be chosen. 

This crystal was silvered uniformly 

on both sides and washed with several organic solvents before mounting in 

its holder. The crystal was placed in a dewar flask and the dewar sealed in 
position in the vacuum system. 


Since the LC circuit across the crystal was a very small value of inductance 
it is reasonable to assume that the value of the voltage induced in it by the 
signal generator would remain constant for the frequency width corresponding 
to the crevasse and experiments confirmed this within !/, percent. 


The experiment consists of tuning the LC circuit to the crystal frequency 
and adjusting the input voltage to a constant value for each of the gases 
studied. A heterodyne frequency meter indicates the frequency of the input 
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signal and the potential across the crystal is read off the galvanometer 
seale. An excellent slow motion drive connectad to the main tuning shaft 
of the signal generator made it possible to change the input frequency very 
gradually. 

The procedure adopted in those measurements was the one suggested by 
Dye. A complete determination of the shape of the crevasse was made from 
observations of the voltage across the crystal at various frequencies on either 
side of the resonance frequency. 

The crystal used in this experiment had a natural fundamental frequency 
of 414,216 ke/see., where a very prominent crevasse was present and a smaller 
creavasse at a frequency about 200 cycles off. Calculations were made on 
the principal dip. 

The main factors contributing to errors were calculated to be of the order 
of + 2,5%. This included fluctuations of the voltmeter, slight creep in gal- 
vanometer and small changes in wave meter frequency, during the course of 
the experiments due mainly to variations in temperature. 

The gases studied here are air, hydrogen, helium, nitrogen, oxygen, argon, 
nitrous oxide and carbon dioxide. All gases with the exception of helium were 
obtained from cylinders and their purity was rated as 99% or better, while 
the helium was obtained from this laboratory was 99,99% pure. 


3. Results 


A typical crevasse for argon is shown in fig. 4 and the values of voltmeter 
and wavemeter readings for all the gases studied are collected in Table I. 
These readings are the mean of several sets of observations taken on different 
days, special precautions being taken to note 20 


the minimum value of voltage, at the bottom i 
of the crevasse. A fluctuation of 2 mm on the “4f~ wu 
galvanometer scale at the bottom of the cre- 16 f 
vasse represents a change of 005 volts only tat E 
across the crystal. In all these experiments, 
the fluctuation of galvanometer deflection at & hay 
the boitom of the crevasse was less than 3mm _ 
constituting an error which was negligible. > at 

The calculations of the equivalent resi- ‚st 
stance, capacitance and Q, the quality factor 
of the crystal were based on the method “7 
suggested by Watanabe®). From the values “Zr ARGON 
of E, the voltage across the crystal at several 
frequencies on either side of the crevasse, one Br ap; 4 CIS KC 
can plot the inverted resonance curve and FREQ. KC/SECOND 
hence, calculate all the electrical constants Fig. 4 


of the crystal including its quality factor Q. 


The calculated values of R’ and K’ together with other relevant physical 
properties of the gases studied are collected in Table II. 


6) Y. Watanabe, Proc. I. R. E. 18, 695 (1930). 
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The shape of the crevasse obtained in these investigations coincides with 
that of Dye. It is seen that the gradient of the voltage across the crystal is 
gradual on the low frequency- side and steep on the high frequency side. The 
frequency width of the crevasse is about one kilocycle while Hubbard’s’) 
observations show a very sharp crevasse. Readjustments of the LC circuit did 
not produce any appreciable changes in either the shape or frequency width 
of the crevasse. 


4. Diseussion 


In ultrasonic interferometry, the shape and peak structure of the crevasse 
is utilised for the determination of sound velocity, absorption and reflection 
coefficients. In these experiments the effect of the medium on the shape and 
equivalent electrical constants of the crystal itself was studied. It is found that 
the medium in which the crystal vibrates does not alter the shape of the 
crevasse. 


Table I 


Vacuum | Air | Hydrogen 
(Temp. 24,2° C) | (Temp. 24,2° C) |. (Temp. 24,0° C) 


Freq. ke/sec | E volts | Freq. ke/sec Evolts _ Freq. ke/sec E volts 


413,370 | | 413,370 413,370 1,77 
413,650 413,425 413.442 1.79 
413,865 413.550 413,532 
413.915 413,618 413,690 
414,012 413,732 413,854 
414,055 413,794 413,914 © 
414.110 413.854 413,956 
414.155 413.907 414.012 
414.175 413.920 414.048 
414.205 413.970 414,090 
414.215 414.000 414.115 
414.221 414.036 414.145 
414.234 414,072 414.163 
414.262 414.108 414.180 
414.275 414.146 414.200 
414.293 414.187 414.210 
414.318 414.204 414.216 
414.342 414.216 414.240 
414.367 414.240 414.253 
414.400 414.264 414.265 
414,451 414,300 414,277 
414.518 414.318 414.307 
414.620 414.336 414.330 
414,747 414.355 414.373 
414.385 414.415 
414,415 414,433 
414.480 414.480 
414.560 414.520 
414.783 414.580 
414.607 
| 414,790 
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?) J. C. Hubbard, Physic. Rev. 41, 523 (1932). 
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Table I contd. 


Carbon Dioxide | Oxygen | Helium 
Temp.: 24,2° C Temp. 24,2° C Temp. 24,0° C 
Freq. ke/sec E volts Freq. ke/sec E volts Freq. ke/sec| E volts 


413,370 1,77 413,370 1,79 413,370 
413,430 413,491 1,82 | 413,430 
413,503 413,673 | 413,476 
413,618 414,794 413,491 
413,645 414,896 413,510 
413,835 414,934 413,562 
413,907 414,012 413,612 
413,975 414,055 413,721 
414,030 414,096 413,769 
414,078 414,120 413,890 
414,102 414,150 413,908 
414,127 414,174 413,932 
414,156 414,187 413,958 
414,169 414,216 413,981 
414,200 414,222 414,025 
414,210 414,246 414,054 
414,216 414,270 414,072 
414,222 414,294 414,102 
414,234 414,312 414,127 
414,258 414,342 414,145 
414,277 414,373 414,156 
414,300 414,403 414,174 
414,325 414,433 414,193 
414,336 414,500 414,205 
414,360 414,536 414,216 
414,400 414,645 414,227 
414,426 414,246 
414,457 414,270 
414,518 414,286 
414,579 414,306 
414,645 414,318 
414,890 414,360 
414,400 
414,408 
414,438 
414,494 
414,548 
414,612 
414,740 
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a) Relation between equivalent electrical resistance and acoustic impedance 

If the effect of the coupled fluid column on the oscillating crystal is treated 
as a parameter of the electrical equivalent circuit, from equations 2 and 3, 
one can see that the value of the equivalent resistance and reciprocal capaci- 
tance should vary linearly with P C provided S and T occurring therein are 
made equal to unity. This condition was secured in these experiments by 
making sure that the sound waves from the crystal surface were completely 
absorbed at the farther end of the container by suitably placed absorbents. 
In fig. 5, the values of P C of the media are plotted against the experimental 
values of R’ the equivalent resistance and it is seen that the relation is linear 
with the range of experimental error. Air and Nitrogen seem to depart from 
this linear relation by about 6%. Probably Nitrogen is the only exception 
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Table I contd. 
Nitrous Oxide Nitrogen Argon 
Temp.: 24,0° C Temp.: 24,0° C | Temp.: 24,8°C 
Freq. ke/sec| Evolts  Freg.ke/see| E volts | Freq.ke/sec | Evolts 
| | 
413,370 1,79 413,370 1,79 | 413,370 1,77 
413,400 1,80 | 413,430 a 413,480 1,76 
413,455 | . 1,82 | 413,491 1,810 413,660 1,73 
413,491 1,82 | 413,562 1,810 414,794 1,65 
413,550 1,82 413,62 | 1,82 413,854 1,62 
413,635 183 | 413,780 | 1,78 413,890 1,55 
413,822 1,78 © 413,854 1,72 413,925 1,32 
413,896 145 413,914 1,00 413,956 0,91 
413,932 126 | 413,932 1,51 413,975 1,34 
413,958 167 | 413,975 1,69 414,007 1,53 
413,975 1,74 | 413,986 1,69 414,067 1,48 
414,000 1,73 | 414,041 1,65 414,127 1,25 
414,048 1,65 | 414,060 1,60 414,145 1,11 
414,097 1,50 414,108 1,26 414,162 0,96 
414,140 1,26 414,140 | 1,11 414,193 0,67 
414,162 1,04 414,162 | 0,87 414,205 0,47 
414,175 0,91 414,195 | 0,59 414,216 0,19 
414,193 0,73 414,200 0,47 414,227 ar 
414,204 0,57 414,210 ae 414,240 0,31 
414,210 0,39 414,216 0,16 414,264 0,62 
414,216 0,18 | 414,227 sit 414,293 0,91 
414,220 0,20 | 414,246 | 0,28 414,325 1,22 
414,228 wt | 414,270 | 0,57 414,355 1,51 
414,258 0,37 | 414,293 | 0,88 414,402 1,63 
414,270 059 | 414,318 1,19 414,457 1,71 
414,290 0,85 | 414,348 | 1,40 414,518 1,74 
414,306 1,00 414,33 | 1,51 414,668 1,71 
414,324 1,19 414,408 1,62 
414,343 1,33 414,435 1,69 
414,376 1,46 414,518 | 1,76 
414,385 1,56 414,536 1,77 
414,404 1,61 414,645 1,78 
414,439 1,69 414,746 1,78 | 
414,463 1,73 414,828 in | 
414,500 1,76 | | 
414,535 1,787 
414,628 1,82 | 
414.205 1.81 | | 
414,820 1,78 | 
Table II 
| Viseosity 
PC|ı K' - 
Medium on | 1/K 24° C Q 
EE es 3. | 4. 5 6. 
Vacuum... .. | 455,0 3,07-10- | 0,325 . 101 27540 
Hydrogen 10,93 512,4 | 2,74-10-™ | 0,364-10% | ,000097 27328 
Helium ..... 16,61 526,4 2,71-10- | 0,369-10% | ,000199 26978 
Nitrogen 40,08 | 582,2 2,56-10-4 | 0,390-10% | ,000186 25785 
| 40,97, 591,8 2,75-10-™ | 0,363-10% | ,000180 23589 
Oxygen ..... | 43.43 655.1. 2.67-10-“ | 0,373-10% | ‚000215 21920 
‚54,56 705,3 | 2,37-10-™ | 0,421-10% | ‚000230 22972 
Carbon dioxide 48,84 649,8 2,31-10-4 | 0,434-10% | ‚000150 25664 
Nitrous Oxide . . 49,67) 649,8 2,54-10-1  0,393-10" | ‚000160 23277 
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and air being a mixture of nitrogen (80%) and oxygen varies as much as 
Nitrogen. 


i a. | 


—R 


Fig. 5 


b) Effect of the acoustic impedance on equivalent capacitance 

The reciprocal of the equivalent capacitance of the quartz crystal at re- 
sonance also should vary linearly with acoustic impedance (equation 5). From 
column (1) and column (4) of Table II ist is found that the variation is not 
linear. 

€) Change of Q with viscosity 

It has been shown®) that in the case of a quartz crystal oscillating in li- 

quids, the quality factor Q decreases exponentially as the viscosity of the 


27+ 
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00006 00008 00010 00012 0001 00016 00078 00020 00022 „00024 
Viscosity poise 
Fig. 6 
8) S. Parthasarathy and A. F. Chhapgar, Ann. Physik 12, 316 (1953). 
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liquid increases. In other words, it means that the damping of the ampli- 
tude of oscillation is more in liquids of greater viscosity. In fig. 6, the Q of the 
crystal used in these experiments is plotted against the viscosity of the media. 
A straight line drawn from Hydrogen, a gas of the least viscosity, to Argon, 
of maximum viscosity (amongst gases barring Neon not studied here), lies 
more or less evenly between the several observed points, the vertical line passing 
through each point representing the limits of experimental error. The only 
exception here is Helium which departs by about 9% from this line. It is 
possible that even in the case of gases the relation between Q and viscosity 
is exponential, but due to the short range of viscosities obtaining in gases, the 
small part of an exponential curve appears linear. 


5. Conelusion 


We conclude from these experiments that the equivalent electrical resi- 
stance of a quartz crystal is proportional to the acoustic impedance of the gas 
in which it is vibrating and that the quality factor decreases exponentially 
with increasing viscosity as in liquids but the range being smaller, the relation 
is linear. 


New Delhi 12, National Physical Laboratory of India. 


Bei der Redaktion eingegangen am 24. Februar 1954. 
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Absorptions- und Leitfähigkeitsmessungen an CdS 
im weichen Röntgengebiet 


Von E. Schnürer 
Mit 6 Abbildungen 


Inhaltsübersicht 


Es wird gezeigt, daß sich die Absorptionskanten der Komponenten des 
CdS in der photoelektrischen Leitfähigkeit dieses Stoffes ausprägen. 

Ziel der vorliegenden Arbeit war, Absorptions- und Leitfähigkeitsmessungen 
im weichen Röntgengebiet durchzuführen und festzustellen, ob sich an den 
Cadmium-K-, Cadmium-L- und an der Schwefel-K-Absorptionskante Un- 
regelmäßigkeiten feststellen lassen. Mit einer Versuchsapparatur, bei der sich 
die zu untersuchenden Kristalle im Röntgenröhrenvakuum befanden, wurde 
die Absorption von Cadmiumsulfid-Einkristallen in der Umgebung der Cad- 
mium-K-Absorptionskante gemessen ; bei den Leitfähigkeitsmessungen wurden 
die theoretisch zu erwartenden Sprünge an den Cadmium-K- und Cadmium-L- 
Absorptionskanten beobachtet. In einer nachträglich bekanntgewordenen 
Arbeit von S. W. Swetschnikow!) wird ein Sprung der Empfindlichkeit 
(Leitfähigkeit) bei A = 0,46 A (Cadmium-K-Absorptionskante) mitgeteilt. Diese 
Arbeit wurde allerdings ohne Filterung der Röntgenstrahlen durchgeführt. 


Versuchsapparatur 


Soweit bekannt, sind bisher im weichen Röntgengebiet an CdS noch keine 
Messungen durchgeführt worden. Um die Absorption der Luft auszuschalten, 
wurde eine Röntgenröhre mit einem seitlichen Ansatz (Meßkammer) gebaut, 
in dem die zu untersuchenden Kristalle angeordnet wurden. Die MeBapparatur © 
ist schematisch in Abb. 1 dargestellt. Da sich zwischen Röntgenröhre und 
Absorber keine Trennwand befindet, wurden die mit den Röntgenstrahlen 
austretenden reflektierten Elektronen mit Hilfe eines starken Magnetfeldes 
am Eindringen in das Meßgefäß gehindert. Die Röntgenstrahlen wurden durch 
entsprechend ausgewählte Kupfer- und Aluminium-Filter homogenisiert. 

Eine Ionisationskammer bzw. ein Interferenzzählrohr, die im Röntgen- 
gebiet zwischen 2 und 5 A brauchbar sind, stand nicht zur Verfügung. 

Deshalb wurde als Indikator eine CdS-Meßzelle am äußeren Ende der Meß- 
kammer angebracht, vgl. Abb. 1. Der Photostrom dieser Zelle dient bei den 
Absorptionsmessungen als relatives Maß für die Intensität der Röntgen- 
strahlen. Für die Leitfähigkeitsmessungen wurde diese Zelle direkt benutzt. 


1) 8. W. Swetschnikow, Techn. Phys. USSR 1952, 22; vgl. auch: Trohn- 
meyer, Glocker, Messner, Naturwissensch. 1958, 12. 
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Der Querschnitt des Röntgenstrahlenbündels wurde so bemessen, daß der 
gesamte Kristall der CdS-Zelle ausgeleuchtet war. Zwischen das Filter, das 
ohne Öffnen der Röhre ausgewechselt werden konnte, und der Indikatorzelle 
wurde der zu untersuchende Absorber gebracht; die eventuell aus dem Ab- 
sorber ausgelösten äußeren Photoelektronen wurden durch ein elektrisches 
Feld abgesaugt. Mit der relativ kleinen Röntgenstrahlintensität konnten nur 
kleine Leitfähigkeitsströme erzielt werden. Aus diesem Grunde wurden sie 
mit einem Gleichstrommeß- 
verstärker gemessen, dessen 
Eingangswiderstand etwa 
10° Ohm und dessen größte 
Empfindlichkeit 10-12 Amp. 
pro Skalenteil betrug. Da 
der Eingangswiderstand des 
Meßverstärkers und die 
CdS-Meßzelle hintereinan- 
der geschaltet waren, wurde 
der Spannungsabfall an der 
Zelle laufend kontrolliert 
und wenn nötig nachge- 
regelt, so daß an der CdS- 
Abb. 1. Schematische Ansicht der Versuchsappa- Zelle immer dieselbe Saug- 
: mmer, B: eld zur ung vaga- . . 
sorberträger und Absorber, 5: Meßzelle, MV: MeBver- in Röntgenmonochroma- 
stärker, i,: Photostrom der Meßzelle, U,: Spannung tor nicht zur Verfügung. 
an der Mebzelle, % : Strom der aus dem Absorber Es wurde zur Homogeni- 
ausgelösten äußeren Photoelektronen sierung der Strahlung ein 
Filterverfahren angewandt, 
so daß dıe Halbwertsbreite nicht mehr als 20%, des kurzwelligen Endes des 
Röntgenbremsspektrums betrug. Z. B. lag bei 
der Röhrenspannung von 33,5 kV 
die kurzwellige Grenze bei 0,37 Ä 
und mit Benutzung eines Filters aus Al 0,lmm Cu 0,1 mm 
war die Halbwertsbreite des homogenisierten R-Bündels . . . 0,3Ä 0,08 A 

Jedes Filter ließ sich in einem größeren zusammenhängenden Wellen- 
längenbereich benutzen, ohne daß sich die Halbwertsbreite wesentlich än- 
derte. Daher ließ sich die Methode, die Messungen in Abhängigkeit von der 
Grenzwellenlänge des Bremsspektrums durchzuführen, rechtfertigen. 

Für die Absorptionsmessungen genügten relative Intensitätsmessungen 
wie sie mit einer CdS-Meßzelle durchgeführt werden konnten. Die in unserem 
Institut hergestellten CdS-Zellen zeichnen sich durch eine außerordentlich 
gute Konstanz aus. 

Messungen an CdS mit einem CdS-Indikator vorzunehmen, wird im ersten 
Augenblick Kritik herausfordern. Die Meßmethode und ihre Anwendung für 
Absorptionsmessungen sind einwandfrei, wenn sichergestellt ist, daß der 
Photostrom einer solchen CdS-Zelle im jeweiligen Wellenlängenbereich 
proportional der Strahlungsintensität ist. Diese Tatsache war auf Grund 
zahlreicher im Institut für Festkörperforschung vorgenommenen Messungen 
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sichergestellt. Darüber hinaus sei 
bemerkt, daß die Kristalle einen 
sehr kleinen Dunkelstrom von 

5 x 10-12 Amp. bei 2 Volt Saug- 

spannung hatten, völlig frei von 
Gleichrichtereffekten waren und 
sich momentan auf den End- 
wert einstellten. 


pi cm") 


x 
Durchführung der Messungen ya 
und ihre Ergebnisse 


a) Die Absorptionsmessungen 


Zwischen 0,3 und 1,5 A wurde ALA} 
von CdS-Einkristallen gemessen. 05 / 15 
In diesem Gebiet war der Photo- Abb. 2. Absorptionsverlauf von (CdS-Ein- 

kristallen zwischen 0,3 und 1,4Ä. (Der Ab- 

sorptionssprung entspricht der Cd-K-Kante.) 

Wertetabelle zur Abb. 2 u: Absorptionskoeffizient, A: Wellenlänge, x: 

gemessene Werte; die ausgezogene Kurve ist der 
berechnete Absorptionsverlauf 


strom der Indikatorzelle 
proportional der Intensität 
der gefilterten Réntgen- 
strahlung. Daher läßt sich 
bei bekannter Dicke des 
Absorbers der Absorptions- 
koeffizient aus relativen 
Messungen bestimmen. In 
Abb. 2 sind die erhaltenen 
Meßpunkte eingetragen. 
Die ausgezogene Kurve gibt 
den berechneten Absorp- 
tionsverlauf in diesem Ge- el 
biet wieder. Mit Hilfe der 6! 1 Pig 
. N 
absorptionskurve wurden 


die Absorptionskoeffizien- Abb. 3. Bestimmung der Lage einer Absorptionskante 

ten vom Cd berechnet, die mit ri _ mem Dem Knick im Photo- 

_  Stromverlauf der e entspricht eine Absorptions- 

ior S wurden aus Lan- 7,nte. Bei Cu liegt die K-Kante vor und bei CdS 

oldt-Börnstein ent die Cd-K-Kante. ip: Photostrom der Meßzelle, U,: 
nommen; daraus erhält Röntgenröhrenspannung, A: Grenzwellenlänge 
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man die Absorptionskoeffizienten von CdS mit Hilfe der Formel: 


o/cas 100\@/s 100 \ /ca 


(&,ca: gewichtsprozentualer Anteil von S und Cd an CdS, «/e Massenabsorptionskoeffi- 


zient.) 

Die für die Absorptionsmessungen ausgesuchten Kristalle sahen homogen 
aus und hatten eine gleichmäßige Dicke. Das Röntgenstrahlenbündel, das den 
Absorber durchsetzte, wurde wie üblich bei Absorptionsmessungen möglichst 
schmal ausgeblendet. 

Die Messungen zeigen deutlich das Auftreten der Cd-K-Kante im Ab- 
sorptionsverlauf der CdS-Kristalle. Die gemessenen Werte liegen durchweg 
höher als die berechneten ; dies wird durch die merkliche Halbwertsbreite der 
Röntgenstrahlung verursacht, die längerwellige Anteile enthält. 


In Abb. 3 wird der Photostrom der Indikatorzelle in Abhängigkeit von 
der Wellenlänge der eingestrahlten Röntgengrenzstrahlung angegeben, wobei 
ein CdS- oder ein Cu-Absorber in den Strahlengang gebracht wurde. Das 
Réntgenlicht wird im Fall von CdS mit 0,1 mm Cu und im Fall von Cu mit 
0,1 mm Al gefiltert, die Knickstellen im Photostromverlauf der Indikatorzelle 
zeigen hier unmittelbar die Lage einer Absorptionskante an. In beiden Fällen 
handelt es sich um K-Kanten (Cd und Cu). 


b) Leitfähigkeitsmessungen 


Mit der gefilterten Röntgengrenzwellenstrahlung wird eine CdS-Zelle ange- 
regt. In einem Bereich, in dem die Halbwertsbreite des homogenisierten 
Strahles 20% des Grenzwellen- 
Ip [Amp] wertes nicht übersteigt, wird das 
| gleiche Filter beibehalten. Die 

Röhrenspannung (und damit die 

Grenzwellenlänge) wird variiert 

und in Abhängigkeit davon der 

Photostrom der CdS-Zelle ge- 

messen. Im Gebiet von 2 bis 7 kV 

Röhrenspannung wird ein Al-Fil- 

ter der Dicke 0,1 mm benutzt und 

im Gebiet von 26 bis 31 kV ein 

Cu-Filter von 0,1 mm. Die 

zur Röhrenspannung gehörenden 

2 28 Up LKV) Grenzwellenlangen sind mit ein- 

+ Pe + getragen. Die beiden folgenden 

an 16 TA) Abbildungen zeigen das Ergehnis 
Abb. 4. Photostromsprung einer CdS-Zelle an der Messungen. 

der Cd-K-Kante. i,: Photostrom, U,: Röntgen- Wir finden in Abb. 4, daß 

röhrenspannung, 1: Grenswellenlinge sich der Photostrom im Gebiet 

von 0,46A bis 0,43A und in 

Abb. 5 im Gebiet von 3,5 A bsi 2,8 A unstetig andert. Er springt beim Uber- 

streichen dieser Gebiete, mit abnehmender Wellenlänge der Grenzstrahlung 

plötzlich auf einen größeren Wert. Bei 3,078 A , 3,319 A und 3,496 Ä liegen die 
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drei L-Absorptionskanten und bei 0,463 Ä die K-Kante des Cd. Da die 
beiden registrierten Unstetigkeiten des Photostromes mit: den L-Kanten 
bzw. der K-Kante des Cd zusammenfallen, wirkt sich hier offenbar der Ab- 
sorptionssprung des Cd auf den im CdS angeregten Photostrom aus. Absorp- 
tionssprung und Photostromsprung verlaufen gleichsinnig. Diese Photostrom- 
sprünge treten am deutlichsten hervor, wenn die CdS-Kristalle möglichst 
dünn sind. 

Im weichen Röntgengebiet wurde die Abhängigkeit des Photostromes der 
CdS-Zellen von der Röntgenstrahlenintensität gemessen. Parameter ist die 
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18 KV (~65A) 
|~ 


‘ 
Up (KV) 0 —i, [mA] 


fi 
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seas Abb. 6. Photostrom einer CdS-Zelle 
Abb. 5. Photostromsprung einer CdS-Zelle in Abhängigkeit von der Röntgen- 
an den Cd-L-Kanten. i,: Photostrom, strahlenintensität (relativ). ir: 
U,: Röntgenröhrenspannung, A: Grenz- Photostrom, iz: Röntgenröhren- 
wellenlänge strom 


Röntgenröhrenspannung. Von 4,5 kV Grenzstrahlung aufwärts wurde stets 
Proportionalität zwischen Photostrom und Röntgenintensität gefunden. 
Zwischen 1,8 kV (6,5 A) und 4,5 kV (2,8 A) Grenzstrahlung waren Photo- 
strom und Röntgenintensität einander nicht proportional. In Abb. 6 ist der 
Photostrom über dem Röhrenstrom (dem die Röntgenintensität proportional 
ist) doppelt logarithmisch aufgetragen. 

Bei 6,5 A Grenzwellenlänge beträgt die Steigung der Geraden !/,; bei 5 Ä 
Grenzwellenlänge beträgt die Steigung erst ?/,, dann folgt ein Knick im Photo- 
strom und bei weiterer Intensitätssteigerung hat die Steigung einen Wert 
größer als 1. Es ist möglich, daß dieses Intensitätsverhalten des Photo- 
stromes auf das Mitwirken der Schwefel-K-Kante. die bei 5 Ä liegt, zurück- 
zuführen ist. 
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Die Ergebnisse der Messungen sind so zu erklären: 

‘ Dort, wo die Absorption eine größere ist, werden aus den entsprechenden 
Termen des CdS mehr Photoelektronen abgelöst und infolgedessen können 
mehr Elektronen in das Leitfähigkeitsband des CdS gelangen; damit wird 
der Photostrom auch größer. Ändert sich nun der Absorptionskoeffizient 
unstetig, wie es an den Absorptionskanten der Fall ist, dann ist zu erwarten, 
daß sich der Photostrom unstetig ändert, so wie es auch gemessen wurde. An 
einer Absorptionskante haben die eingestrahlten Röntgenquanten gerade die 
Energie, die einem bestimmten Elektronenniveau entspricht. Das heißt, 
die Energie der Röntgenquanten reicht gerade dazu aus, um ein Elektron von 
einem solchen Term abzulösen und es in das Leitfähigkeitsband zu befördern. 
Ist die Energie der Röntgenquanten größer als die einem solchen Term ent- 
sprechende, dann bekommt das abgelöste Elektron die Überschußenergie als 
kinetische Energie mit. Diese Überschußenergie kann entweder an das Gitter 
des festen Körpers abgegeben werden, oder Sekundäreffekte verursachen; 
ein Elektron mit genügend hoher kinetischer Energie kann auch den Verband 
des festen Körpers als äußeres Photoelektron verlassen. (Mit den bisherigen 
Methoden konnten noch keine äußeren Photoelektronen nachgewiesen werden.) 
Im Bereich geringer Absorption bestimmen diese Sekundäreffekte (ausge- 
löste Sekundärelektronen, sekundäre Röntgenstrahlung, Anregung von Rönt- 
genfluoreszenzstrahlung usw.) vorwiegend die Leitfähigkeit des CdS. Die 
Wahrscheinlichkeit für das Auftreten dieser Sekundäreffekte wird mit zu- 
nehmender kinetischer Energie der abgelösten primären Photoelektronen 
sehr groß. 

Da die Absorption von CdS bei Röntgenanregung vorwiegend durch die 
Absorption des Cd bestimmt wird, stammen die beim Absorptionsakt primär 
abgelösten Photoelektronen vorwiegend aus Cd-Termen. Lediglich im Bereich 
der Cd-L-Kanten und im Bereich der S-K-Kante (des großen Absorptions- 
sprunges wegen) ist die S-Absorption für die Gesamtabsorption des CdS mit- 
bestimmend. Bei der Rekombination werden Löcher in den tieferen Termen 
von Elektronen aus darüber liegenden Termen wieder aufgefüllt (Lebens- 
dauer von Löchern in tieferen Termen etwa 107 sec. — siehe Wien-Harms, 
Handb. d. Exp. Phys., Bd. 24, 1); bei dieser Rekombination wird die ent- 
sprechende Röntgenfluoreszenzstrahlung emittiert. Sind dann die Löcher 
im Valenzband angelangt, dann setzt der Rekombinationsmechanismus ein, 
wie man ihn von der Anregung mit sichtbarem Licht her kennt. Löcher in 
tiefen Termen können wegen ihrer kurzen Lebensdauer und infolge ihrer Unbe- 
weglichkeit keinen Beitrag zum Photostrom liefern. 

Die vorliegenden Messungen haben zunächst nur orientierenden Wert. 
Die Untersuchungen werden nach Fertigstellung eines Röntgenmonochro- 
mators fortgesetzt. 


Herrn Prof. Dr. F. Möglich und Herrn Dr. H. Simon danke ich für ihre 
Unterstützung und für anregende Diskussionen. 


Berlin-Buch, Institut für Festkörperforschung der Deutschen Akademie 
der Wissenschaften. 


Bei der Redaktion eingegangen am 6. März 1954. 
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Zur Berechnung von Elektroneninterferenzen 
an Zwillingskristallen’) 


Von Heribert Jahrreiß 
Mit 4 Abbildungen 


Inhaltsübersicht 


Es wird untersucht, inwieweit die elektronenoptischen Phasenbeziehungen 
zwischen den Komponenten eines Zwillingskristalls für die Berechnung des 
Beugungsdiagramms eine Rolle spielen. An einem speziellen Beispiel (flächen- 
zentriert-kubisches Gitter, verzwillingt nach (111)) wird gezeigt, daß diese 
Phasenbeziehungen vernachlässigt werden können, falls der Kristall in Rich- 
tung senkrecht zur Zwillingsebene die Dicke weniger Elementarzellen — etwa 
30 bis 40 A — überschreitet. Methode und Ergebnisse sind auf andere Gitter- 
konfigurationen übertragbar. 


In einer früheren Arbeit?) wurde bereits das Problem der elektronen- 
optischen Phasenbeziehungen zwischen den Komponenten eines Zwillings- 
kristalls angedeutet: Man berechnet das Elektroneninterferenzdiagramm 
eines aus mehreren Kristalliten bestehenden Präparates üblicherweise so, daß 
man die Intensitätsverteilungen (d. h. die Quadrate der Beugungsamplituden) 
der einzelnen Kristallite überlagert oder mit anderen Worten die Ewaldsche 
Konstruktion anwendet auf ein reziprokes Gitter (rezG), das eine Ineinander- 
schachtelung der rezG der einzelnen Kristallite darstellt. Die Berechnung 
der Intensitätsverteilyrrgen ist dabei in bekannter Weise unter Berücksich- 
tigung der üblichen Faktoren durchzuführen. Dieses Vorgehen ist streng 
richtig, wenn die einzelnen“Kristallite unabhängig voneinander streuen, also 
keine elektronenoptische Phasenbeziehung zwischen ihnen besteht. Wir 
wollen diesen Fall den ‚„Inkohärenzfall‘‘ nennen. — Bestehen dagegen zwischen 
einzelnen Kristalliten solche Phasenbeziehungen, etwa dadurch, daß zwei 
räumlich verschieden orientierte Kristallite der gleichen Substanz eine Netz- 
ebene gemeinsam haben (Verzwillingung), so sprechen wir vom „Kohärenz- 
fall“. Hier sind die abgebeugten Amplituden der beteiligten Kristallite zu 
überlagern und aus deren Summe durch Quadrierung die Gesamtintensitäts- 
verteilung zu bestimmen. Die Ineinanderschachtelung der rezG der Einzel- 
komponenten gibt also die Verhältnisse im Kohärenzfall nicht genau, höch- 


1) Vorgetragen auf der Frühjahrstagung 1954 der Nordwestdeutschen Physikalischen 
Gesellschaft in Goslar. 
2) H. Jahrreiß, Ann. Physik (6) 14, 319 (1954). 
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stens näherungsweise richtig wieder. Da aber diese Methode bisher ver 
schiedentlich bei der Diskussion von Beugungseffekten an Zwillingskristallen 
benutzt wurde (Cochrane°), Kirehner und Rüdiger®), Menzer®), vgl. 
auch den zusammenfassenden Bericht von Raether®), S. 91ff.), soll die 
Berechtigung dieses Verfahrens im folgenden gepriift werden. 

Dazu muß man das Ergebnis von Berechnungen der Intensitätsverteilung 
der gleichen Kristallkonfiguration mit und ohne Berücksichtigung der Phasen- 
beziehung miteinander vergleichen. Den Berechnungen wird ein Modell 
zugrunde gelegt, dem einige idealisierende Bedingungen vorgeschrieben 
werden müssen. Trotz dieser Einschränkung sollte das Ergebnis auch die 
Beurteilung der Verhältnisse im allgemeinen Fall mit hinreichender Genauigkeit 
gestatten. — Es wird sich zeigen, daß ein nennenswerter Unterschied zwischen 
Kohärenz- und Inkohärenzfall nur auftritt, wenn beide Kristallite in der 
Richtung senkrecht zu der gemeinsamen Netzebene nur wenige Elementar- 
zellen übereinander aufweisen, d. h., wenn wir es näherungsweise mit Flächen- 
gittern zu tun haben. Daher reicht die geometrische Theorie der Elektronen- 
interferenzen als Rahmen der folgenden Überlegungen aus. 


IL. 


Die weiteren Betrachtungen beziehen sich auf ein spezielles Modell mit 
folgenden Eigenschaften: Man habe zwei miteinander verzwillingte Gitter 
beliebiger Art. Diese sollen senkrecht zur Zwillingsebene gleichviel Atomlagen 
übereinander aufweisen, wobei je eine zur Zwillingsebene parallele Ebene die 


obere und untere Begrenzung des gesamten M Atomlagen aufweisenden 
Komplexes bilde. Nach den beiden anderen Raumkoordinatenrichtungen 
sei der Komplex beliebig weit ausgedehnt. Die beiden Zwillingskomponenten 
seien mit I und II bezeichnet. — Es ist bequem, für die weiteren Überlegungen 
zu einem Koordinatensystem überzugehen, dessen eine Achse die Richtung 
senkrecht zur Zwillingsebene habe (,‚c-Achse‘‘); die beiden anderen Rich- 
tungen weisen dann in der Zwillingsebene von einem Atom zu je einem seiner 
nächstbenachbarten Atome. In der c-Richtung entspricht die Identitäts- § 
periode einer bestimmten Anzahl n von Atomlagen übereinander. Innerhalb 
des neuen Koordinatensystems gibt es für I und II je eine Elementarzelle; 
die Schichthöhe von I und II in c-Richtung entspricht dann je M/2n = M 
Elementarzellenhöhen. Im neuen Koordinatensystem erhalten die einzelnen 
Netzebenenscharen natürlich andere Indices. — Wir haben nun für dieses 
Kristallitmodell die Intensitätsverteilung zu berechnen einmal mit und einmal 
ohne Berücksichtigung der Phasenbeziehung: 

1. Für ein allseitig unendlich ausgedehntes Gitter erhält man ein rezG 
mit punktförmigen Intensitätsbereichen. Die Intensität in den intensitäts- 
behafteten Punkten bestimmt der Strukturfaktor. Wir betrachten nun die 


. Gittergeraden des rezG, die senkrecht zur Zwillingsebene verlaufen, d.h. die 


geometrischen Orte der Gitterpunkte des rezG mit den Indices (a, a, =), 
wo a, und a, ganze Zahlen, x laufende Koordinate. Haben I und II — wie 


3) W. Cochrane, Proc. physic. Soc. 48, 723 (1936). 

4) F. Kirchner u. O. Rüdiger, Ann. Physik (5) 30, 609 (1937). 
5) G. Menzer, Z. Kristallogr. 99, 378 (1938) u. 99, 410 (1938). 
6) H, Raether, Ergebnisse d. exakt. Naturw. XXIV, 54 (1951). 
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ursprünglich vorausgesetzt — eine nur endliche Ausdehnung in Richtung 
senkrecht zur Zwillingsebene, während die beliebig weite Ausdehnung in den 
beiden anderen Richtungen erhalten bleibt, so erweitern sich die Intensitäts- 
bereiche längs dieser Gittergeraden, während senkrecht dazu keine Bereich- 
erweiterung stattfindet (Bragg und Kirchner’), sowie Kirchner®)). Die 
Problemstellung reduziert sich also auf die Frage nach der Inten- 
sitätsverteilung längs dieser Geraden. Ist nämlich diese Intensitäts- 
verteilung bekannt, dann ist das resultierende Beugungsdiagramm sofort 
durch Anwendung der Ewaldschen Konstruktion zu ermitteln. Bei diesem 
Übergang vom rezG zum Beugungsdiagramm kommen natürlich nur Gitter- 
punkte bzw. Intensitätsbereiche des rezG in Frage, die so niedrig indiziert 
sind, daß ihnen im Kristallgitter ausgebildete Netzebenenscharen entsprechen. 
Ferner ist für die Beurteilung der Intensität der einzelnen Reflexe der Flächen- 
häufigkeitsfaktor, Lorentzfaktor usw. in üblicher Weise zu berücksichtigen. 
Uns soll aber im weiteren Verlauf der Überlegungen nur der Intensitätsverlauf 
längs der erwähnten Gittergeraden interessieren, wie er sich aus der Berück- 
siehtigung von Struktur- und Gitterfaktor allein ergibt. 

2. Zur Berechnung der Intensitätsverteilung ohne Berücksichtigung 
der Phasenbeziehung (Inkohärenzfall) haben wir für jede der Komponenten | 
und II die aus dem Strukturfaktor errechenbare Intensitätsverteilung S?(x) 
längs der betreffenden (a, a, x)-Gittergeraden zu multiplizieren mit der vom 
Gitterfaktor herrührenden Intensitätsverteilungsfunktion @%(x, M), so daß 
sich für I und II je eine Intensitätsverteilung £7(x, M) = Sf(x) -Gj(x, M) 
bzw. Zfı(x, M) = Sfi(x)-Gfı(z, M) ergibt. Hier wird speziell Gj] = Gj; wegen 
der Voraussetzung gleicher Anzahlen von Elementarzellenschichten über- 
einander in I und Il. Die Gesamtintensitätsverteilung für den Inkohärenzfall 
wird dann 


(x, M) = Zf(x, M) + (x, M) = {S¥ (x) + Str (x)} -G? (x, M). (1) 
Hierin ist zu setzen fiir jede Gittergerade 
G%(x, M) = rn (vgl. z. B. v. Laue®), S. 112). (2) 

Hingegen sind die S? und Sf; verschieden für die einzelnen Gittergeraden. 
Bei jedem Gittertyp gibt es aber immer nur ein endliche Anzahl Arten von 
Gittergeraden des rezG in c-Richtung, die sich in ihrer Intensitätsbelegung 
unterscheiden (beim kubisch-flächenzentrierten z. B. 3 Arten). Beim Fort- 
schreiten innerhalb des rezG in a,- oder a,-Richtung kehren die Vertreter 
der einzelnen Arten in periodischer Reihenfolge wieder. 

3. Die Berechnung der Intensitätsverteilung mit Berücksichtigung der 
Phasenbeziehung (Kohärenzfall) soll in der Weise erfolgen, daß die gesamte 
Kristallitschicht betrachtet wird als Aneinanderreihung von Prismen mit der 
gleichen Grundfläche, wie sie die Elementarzellen von I und II aufwiesen, 


und der Höhe der gesamten Schicht — also M = 2n M Atomlagen über- 
einander — als Prismenhöhe (Abb. 1). Jedes solche Prisma mit den darin 
enthaltenen Atomen kann man betrachten als eine ‚‚Mammutelementarzelle‘‘, 
deren Strukturfaktor sich nach Kenntnis der Atomkoordinaten im Inneren 


1) W. L. Bragg u. F. Kirchner, Nature 127, 738 (1981). 
8) F. Kirchner, Ann. Physik (5) 18, 38 (1932). 
9) M. v. Laue, Materiewellen und ihre Interferenzen, 2. Aufl., Leipzig 1948. 


| 
er ver 
vgl. 
soll die 
hasen- 
Model] 
hrie ben 
ich die 
uigkeit 
vischen 
in der 


24 Annalen der Physik. 6. Folge. Band 15. 1954 


q des Prismas berechnen läßt. Dieser Strukturfaktor S? ist nicht mehr wie 
} i oben S? und Sf; nur eine Funktion von x, sondern auch abhängig von der An- 


zahl der Atomlagen M 4 Da der Kristallit nur aus einer Lage dieser Mammut- 
elementarzellen besteht, wird 
; der Gitterfaktor konstant 
1 G* =1, und der Struktur- 
; i faktor gibt für sich allein 
a i die Intensitätsverteilung im 
Kohärenzfall 
Um unsere Rechnungsergeb- 
nisse vergleichen zu kénnen 
mit den Resultaten von 2., 
beschränken wir uns auf 


i's solche M, die durch 2 n teil- 
.1. Zugrunde liegendes ell. 

: Inkohärenzfall (Summierung über je M Elementar- bar sind, so daß gilt 

; zellen in I und II), rechts Kohärenzfall (,,Mammut- Zi (x, M) = S*%(x, M) 
elementarzelle‘‘ mit M = 2n M Atomlagen) mit M = M/2 n. 


IM. 
i Das skizzierte Programm wurde insbesondere durchgefiihrt fiir zwei miteinander 
nach (111) verzwillingte kubisch-flichenzentrierte Gitter. Hier wird speziell n = 3; 
die einzelnen Atomlagen sind von gleicher Belastung und äquidistant. Im neuen Koor- 
’ dinatensystem vom hexagonalen Typ erhalten die der Zwillingsebene parallelen Oktaeder- 
Br: ebenen die Indices 003, 006, .... Man kann nach ihrer Intensitätsbelegung 3 Arten von 
Gittergeraden unterscheiden; die Geraden mit (00x), (10x) und (01x) gehören zu je einer 
Er: a = dieser Arten. Allgemein gelten für die 
(00x) Zugehörigkeit einer Gittergeraden mit 
den Indices (a, a, x) zu einer der 3 
(10x) ——»— Arten die Regeln: 


(01x) Ist für Ia,—a,=m mod 3 oder 
für II a, — a, = m mod 3, so ge- 

hört die betreffende Gittergerade zu 

Be; einer Art, bei der sich für unendlich 
ausgedehnte Kristallgitter Intensitäte- 
belegung findet für ganzzahlige x mit 

| x=m mod 3. Dabei sind also die 
Gittergeraden des rezG von I mitm =0 
identisch mit denen von IT mit gleicher 
Restklassenzugehörigkeit. Außerdem 
decken sich die Gittergeraden von I mit 
3 m = 1 mitdenen von II mitm = 2 und 
Bi umgekehrt. Die Intensitätsbelegung der 
ig . einzelnen Geradentypen zeigt Abb. 2, 
i die Lage ihrer Fußpunkte in der zur 
| Zwillingsebene parallelen Grundebene 
des rezG veranschaulicht Abb. 3. Man 
sieht, daß eigentlich nur mehr 2 Arten 
von Gittergeraden übrigbleiben nach 


(4) 


| ' 0 1 2 3 4 5 6 7x 
; Abb. 2. Intensitätsbelegung verschiedener Überlagerung der rezG von I und II, 
Gittergeraden des rezG für unendlich ausge- wenn man keinen Wert darauf legt 
dehnte Gitter in Grundorientierung (I) und zu unterscheiden, welche Interferenz- 
illingsorientierung (II), sowie für Über- punkte von I und welche von II stam- 

rung beider Typen (I + LI) men. 
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Im Inkohärenzfall erhält man fiir die Gerade mit (00x): 
= Si,(x) = 3+ 4cos}2ax + (5) 
fir die Gerade mit (10x): 
(6) 

A= 
B= sinaax-costax 
D= cos x+ cosh ax, 


ebene parallelen Grundebene. Die Intensitätsverteilung in Abhängigkeit von x für unend- 

lich ausgedehntes Gitter ist bei allen durch ein Quadrat gekennzeichneten Geraden die- 

selbe wie für die (00x)-Gerade in Abb. 2, entsprechend bedeutet ein Dreieck den Typ 
(10x) und ein Stern den Typ (012) 


und mit den gleichen Bezeichnungen 

+2AB-2CD+2C0-—2D, (7) 
für die Gerade mit (01x): 

(x) von (01x)-Gerade = S},(x) von (10x)-Gerade und umgekehrt. (8) 
Für die Geraden mit (10x) und (01x) wird also mit (6) und (7) für die in (1) auftretende 
Summe S32 (x) + =2(1 + AP + BP +0? + 2D) (9) 

= 4(1 + sin? (4.2 x) — cosa x- cos} 2). 

Man beachte die Periodizität aller S?(x) im Intervall 0 <3=n. 


Im Kohärenzfall ergibt die Durchrechnung für M = 0 mod 3 für die Gerade mit 
(a, a, x), wenn wir die Variable zunächst mit x’ bezeichnen: 


M 
32’ 4 , 
expine 


, 
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Diese Funktion ist periodisch im Intervall 0 < 2’ < M. Diesem Intervall entspricht aber 
auf der Gittergeraden des rezG die gleiche Länge, die oben vom Intervall O0 <x<3 
eingenommen wurde. Denn mit der Vergrößerung der Längeneinheit der c-Achse der 
Elementarzelle im Kristallgitter von 3 auf M Atomlagen verkleinert sich gleichzeitig 
die entsprechende Längeneinheit des rezG auf den 3/Mten Teil. Setzt man also die 
Relation x’ = M x/3 = 2 M x in (10) ein und geht über zu S?, so kommt zunächst für 
das Quadrat des Ausdrucks in der geschweiften Klammer (vgl. z. B. v. Laue®) a. a. 0.); 
M-ı 2 in? 

sin? 2 x 

d.h. dieser Term wird nach (2) gleich dem Gitterfaktor G?(x, M), der im Inkohirenzfall 
auftrat, und es ergibt sich eine formale Parallelität zwischen der Strukturfaktorensumme 
des Inkohärenzfalles einerseits und dem Quadrat des in eckige Klammern gesetzten 
Terms von (10) andererseits, weswegen wir diesen im folgenden ,,Pseudostrukturfaktor* 
P?(x, M) nennen wollen. Dieser Faktor wird also beim Übergang von x’ nach x und unter 
Vertauschung einiger Glieder: 


.(2 1 
+ ga) 


+ exp22i3-3 M (12) 


m+ @M + 1) 


2 
(8M + 2)| 


Nach einigen elementaren, aber etwas umständlichen trigonometrischen Umformungen 
erhält man aus (12) für die einzelnen Gittergeraden des rezG verschiedene Ausdrücke für 
P:(x, M): 


für die Gerade mit (00x): 


= 6 + 2 {4 cos20 + x 


+ cos 22 (3 M — 2) + 2000225 (3 M — 1) 


+ 3000277 3M + (8M + 1) (18) 
+ OM +2), 
fiir die Gerade mit (10 x): 
P2(x,M) = 4 x M + cost xx M cost + Bint M «sin? 
+ 2 M- sin [sin (M + 1) 2—sin (M —1) 22] (14) 
— cos x M cos 2-5 [cos (M —1) 2x + cos (M + 1) aah. 
Ferner ergibt sich aus (12) für die Gerade mit (01x), daß immer gilt 
ER (x, M) = Pre (—x, M). (15) 


Mit Hilfe der gefundenen Formeln für I}; und 2? wurden die Kurven- 
züge im sich periodisch wiederholenden Intervall 0 < x <3 für verschiedene 
Werte von M berechnet. Abb. 4 zeigt jeweils beide Kurven für die Intensitäts- 
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verteilung längs der Geraden mit (10 x) mit den Parameterwerten M = 1 
und M = 6. Man sieht deutlich, daB schon für M = 6 der Unterschied zwischen 
beiden Darstellungen relativ zur Intensität im Maximum kaum noch nennens- 
wert, auf der Photoplatte sicher nicht mehr nachweisbar ist, zumal ja auf 


dem Präparatträger immer 
eine Anzahl von Kristalli- 
ten gleichzeitig vom Elek- 
tronenstrahl erfaßt wird, 
deren Dicke — also deren 
M — innerhalb eines ge- 
wissen Bereiches streut; es 
kommt also zu einer Ver- 
wischung der feinen Unter- 
schiede der Nebenmaxima. 

Bemerkenswerterscheint 
mir, daß die Unterschiede 
nicht über den Punkten der 
Abszissenachse auftreten, 
die den Intensitätsmaximis 

unendlich ausgedehnter 

Gitter entsprechen. 

Diese Regel ist für die Git- 
tergeraden mit (10x) und (01x) 
ohne weiteres erfüllt, nicht je- 
doch für diejenige mit (00x), 
wo die aus den Formeln für 
den Kohärenz- und den Inko- 
härenzfall sich ergebenden 
Werte in x=0, x=3 usw. 
zueinander im Verhältnis 2:1 
stehen. Dies liegt an der Eigen- 
tümlichkeit der Theorie, zwei 
Gittern mit M, und M, Atom- 
lagen Intensitätsmaxima zu- 
zuordnen, deren numerische 
Werte zueinander im Verhält- 
nis Mj: M3 stehen. So kommt 
es, daß im Inkohärenzfall durch 
Addition ein Maximumwert 
von 2(M/2)? entsteht, während 
im Kohärenzfall, wo mit einem 


L? 


0 / 2 jx 
Abb. 4a. 53, (gestrichelt) und 2%, (ausgezogen) für 
Gittergerade des rezG (10x) und M=1 sowie Lj, 
(punktiert) fiir gleiche Gerade und M = 6. Alle Kurven 
normiert auf gleiche Höhe für x = 1 und x = 2, Inten- 
sitätswerte in relativen Einheiten 


0 1 2 ix 
Abb. 4b. 23, (gestrichelt) und 2% (ausgezogen) für 
Gittergerade des rezG (10x) und M = 6, zehnfach 
überhöht gegenüber Abb. 4a 


einheitlichen Kristall von M Atomlagen gerechnet wird, ein Maximumwert von M? 
resultiert. Da es hier um die Relationen zwischen den Intensitätswerten an verschie- 
denen Punkten der Gittergeraden geht, erscheint es gerechtfertigt, im Falle (00x) durch 
geeignete Normierung die Maxima beider Kurven einander anzugleichen. — Die Schwie- 
rigkeit tritt im Falle der beiden anderen Gittergeraden nicht auf, weil hier die Maxima 
von I und II an verschiedenen Punkten liegen. 


IV. 


Wir sahen, daß die Unterschiede zwischen kohärenter und inkohärenter 
Darstellungsart mit wachsendem M immer geringer werden. Diese Erkenntnis 
läßt sich natürlich durch Ausrechnung für eine größere Anzahl von M-Werten 
immer weiter erhärten. Doch ist dieses Verfahren insbesondere für höhere ° 
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Schichtzahlen recht mühsam, da wegen der zunehmenden Oszillation der 
Kurven immer mehr Zwischenpunkte berechnet werden miissen, und es bringt 
keine allgemeingültige quantitative Aussage über die Abhängigkeit der beiden 
Darstellungen von M. Es wurde daher nach einem geeignet definierten Ab- 
weichungsmaß gesucht, das nur noch von M abhängig sein sollte. 

Die Intensitätsverteilungen £,(x, M) und S?(x, M) sind periodische 
Funktionen von x mit dem Periodizitätsintervall 0 < x <n. Wir können 
daher einen Vergleich über den Gesamtverlauf beider Funktionen anstellen, 
wenn wir nur ihren Verlauf innerhalb eines Periodizitätsintervalls verfolgen, 
etwa innerhalb 0 < x <n. Wir bilden zunächst die Differenz 


®(x, M) = Z},(x, M) — Zi (x, M). (16) 

Das Quadrat dieses ®(x, M) können wir über eine Periodenlänge nach x inte- 
grieren, und wir definieren in Anlehnung an die Methode der kleinsten Qua- 
drate als Abweichungsmaß die Funktion 
(x, M) dx 

(x, M) dx (17) 


F(M) (u, M) du, (18) 


wenn wir von x übergehen zur Variablen u = x(2a/n) mit 0< u << 2a. — 
Der Faktor 1/M* dient zur Normierung, da 5? und ?, Maxima vom Wert M: 
aufweisen, © also proportional M4 bzw. M? ist. 


Für den Spezialfall des flächenzentriert-kubischen Gitters erhält man als ®2(u, M) 
nach einigen Umformungen und unter Vernachlässigung konstanter Zahlenfaktoren für 
die einzelnen Gittergeraden des rezG verschiedene Ausdrücke: 


für die Gerade mit (00x): 
®: (u, M )= 

für die Gerade mit (10x): 
®: (u, M) = 


sin’ (u M/2) 
(1+ 8 cos2u- (19) 
sin® (u M/2) 
sin?3u 


sin? u cos? (20) 
\ 6 M)| 

Man beachte, daß diese Funktionen verschwinden für die Maxima wie für 
die Minima der auf beiderlei Art gerechneten Intensitätsverteilungsfunktionen, 
die übrigens auch ihrerseits in der Lage ihrer Nullstellen übereinstimmen. — 
Die bestimmten Integrale, die sich durch Einsetzen dieser ©? in (18) ergeben, 
können mit Hilfe des Residuensatzes gelöst werden (s. z. B. Whittaker u. 
Watson!), S.112f.). Man führt dazu eine neue Variable z = expi x ein, 
wobei der Integrationsweg in den Einheitskreis auf der komplexen Zahlen- 
ebene übergeht. Im Falle des flächenzentriert-kubischen Gitters erhält man 
schließlich für die Gerade mit (10 x): 

% A 
F (M)qo2) 

10) E.T. Whittaker u. G. N. Watson, A Course of Modern Analysis, 4th Ed., 

Cambridge 1950. 
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Diese Abhängigkeit des Abweichungsmaßes von M~* bedeutet praktisch, 
wenn man einen Intensitätsempfindlichkeitsspielraum der Photoplatte von 
zwei Zehnerpotenzen zugrunde legt, daß eine Berücksichtigung der 
Kohärenz nur bei extrem dünnen Schichten (M < 5) notwendig ist. 

Freilich scheinen die Verhältnisse für die Gerade mit (00x) weniger günstig zu liegen. 
Hier tritt neben dem M-*-Glied noch ein M-!-Glied auf. (Ersetzt man Prove) durch 

— M/2) 
(00) = Sint Bu” 

das bis auf den Zahlenfaktor 9? = 81 eine Majorante von ®? darstellt, so erhält man für 
M = 0 mod 3 


(22) 


= 1 1/4 6 

M) du = 55 + ) (23) 
Es ist aber zu berücksichtigen, daß auf der (10x)-Geraden für kleine M Unterschiede in der 
Lage der Maxima auftraten (vgl. Abb. 4, M = 1), was auf der (00x)-Geraden für kein 
M vorkommt, denn (5) und (13) enthalten nur cos-Glieder mit relativ zu den Maxima 
symmetrischem Verlauf. Der Unterschied beider Berechnungsarten für die (00x)-Gerade 
liegt somit vor allem im Auftreten der Mten Oberwelle und einiger Wellen mit benach- 
barten Frequenzen im Pseudostrukturfaktor P? und damit in verschiedener Halbwerts- 
breite der Maxima, deren Lage gegenüber den Verhältnissen beim allseitig unendlich weit 
ausgedehnten Kristall nicht verändert wird. Der langsamere Abfall des Abweichungs- 
maßes mit M- ist daher bei der (00x)-Geraden längst nicht so kritisch, wie er es bei der 
(10x)-Geraden wäre. Man wird also auch hier schon von recht dünnen Schichten an auf- 

wärts auf eine Berücksichtigung der Kohärenz verzichten können. 

An dieser Abschätzungsmethode könnte bedenklich erscheinen, daß ja 
grundsätzlich das Kleinwerden der Fläche zwischen zwei Kurvenzügen nicht 
ausschließt, daß bei sonst weitgehender Annäherung an einigen wenigen 
Stellen starke Abweichungen geringer Halbwertsbreite eintreten können. 
Dagegen spricht hier aber die Beschränktheit und Stetigkeit der trigonome- 
trischen Funktionen, aus denen sich die Zähler der ©? zusammensetzen. Die 
einzigen Stellen, an denen die ©? große Werte annehmen könnten, sind die 
Umgebungen der Nullstellen des Nenners. An diesen Stellen werden aber die 
@ selbst gleich 0. 


V. 


Eine Anwendung der hier skizzierten Methode der „Mammutelementar- 
zelle‘‘ auch auf kompliziertere Gitter dürfte keine grundsätzlichen Schwierig- 
keiten bereiten. Die dann auftretenden Abweichungsfunktionen zwischen den 
Intensitätsverteilungen im Kohärenz- und Inkohärenzfall wird man mittels 
eines der folgenden Integrale nach der Majorantenmethode abschätzen können 
(vgl. auch (23)), soweit nicht eine explicite Ausrechnung in Frage kommt: 


1 sin? N x dx 
Ne — = tra (24) 


sin? 
0 


n 
sin N xdx _2n 


(25) 


an 
(26) 


sin? x 2N3 
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Grundsätzliche Schwierigkeiten bei der Abweichungsabschätzung bereitet 
allerdings die Verallgemeinerung auf den Fall, daß die Komponenten I und II 
verschieden dick sind, weil dann kein einheitlicher Gitterfaktor auftritt. 
Dieser Fall ist aber nur interessant, wenn die Dicke von I und II noch einiger- 
maßen vergleichbar ist, weil andernfalls eine Komponente intensitätsmäßig 
vollkommen überwiegen würde. Sind aber I und Il nicht allzu verschieden, 
so wird die Rechnung mit dem symmetrischen Modell eine gute Näherung 
darstellen. 

Als Ergebnis der Überlegungen läßt sich also feststellen, daß eine Berück- 
sichtigung des Kohärenzverhaltens bei Zwillingsbildung nur dann notwendig ist, 
wenn man es mit extrem dünnen Schichten von der Dicke weniger Elementar- 
zellenhöhen zu tun hat. Oberhalb dieser Dicke von größenordnungsmäßig 
etwa 30 bis 40 Ä darf man also, soweit nicht die dynamische Theorie heran- 
gezogen werden muß, in üblicher Weise auf der Grundlage der Laueschen 
Interferenzbedingungen die Beugungsdiagramme der beteiligten Komponenten 
berechnen und additiv zum Beugungsdiagramm der Zwillinge zusammenfügen. 

Unterhalb dieser Grenze aber ist die auf der Gittervorstellung basierende 
Theorie ohnehin nur mit Vorsicht anzuwenden, wie die Arbeiten von Boersch®) 
König!2), Germer und White?) sowie Berry!#) gezeigt haben, in denen 
die Verfasser die Interferenzdiagramme kleiner und kleinster Kristallite ver- 
schiedener Gestalt auf der Grundlage der Debyeschen Gasinterferenz- 
methode 5) berechnet haben, wobei sich teilweise beachtliche Abweichungen 
gegenüber den Voraussagen der geometrischen Theorie ergaben. 

Im allgemeinen wird man im Realkristall nicht die hier zugrunde gelegte 
Struktur von Lamellen finden, die nach zwei Richtungen hin sehr weit ausge- 
gedehnt sind. Die Berechnungen sind aber auch auf diesen Fall übertragbar, 
da die Intensitätsverteilung für die einzelnen Gittergeraden des rezG immer 
nur von der Ausdehnung des Kristallgitters jeweils in Richtung dieser Gitter- 
geraden abhängt. 


u) H. Boersch, Z. Physik 119, 54 (1942). 

12) H. König, Optik 3, 201 (1948). 

13) H. L. Germer u. A. H. White, Physic. Rev. 60, 447 (1941). 
14) Ch. R. Berry, Physic. Rev. 88, 596 (1952). 

15) P. Debye, Ann. Physik (4) 46, 809 (1915). 


Köln, I. Physikalisches Institut der Universität. 


Bei der Redaktion eingegangen am 25. Mai 1954. 
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Magnetische Pulvermuster auf Martensit 
Von Wilfried Andrä 


Mit 8 Abbildungen 


Die Methode der Bitterschen Streifen ist seit einer Reihe von Jahren zu 
einem der wichtigsten Hilfsmittel geworden bei der Klärung von grundsätz- 


lichen Fragen des Ferromagnetismus!). 


Wir haben bei der Anwendung der Methode auf martensitischem Stahl 
interessante Muster gefunden, die für den Problemkreis der martensitischen 


Umwandlungen von Bedeutung sein 
können. Als Probenmaterial wurde 
Stahl der ungefähren Zusammen- 
setzung 1,4% C; 0,4% Mn; 1,7% Cr; 
0,2% V; 0,3% Si; Rest Fe verwendet, 
in dem sich nach geeigneter Wärme- 
vorbehandlung vereinzelte große Mar- 
tensitnadeln innerhalb von Restau- 
stenit gebildet hatten. Danach wurde 
die Probe angelassen und anschlie- 
Bend durch Tiefkühlung nochmals 
Martensit erzeugt, so daß tetragonaler 
und angelassener Martensit nebenein- 
ander vorhanden waren. Ein Schliff- 
bild des Materials zeigt Abb.1. Die 
Härteeindrücke dienen zur Markie- 
rung. In Abb. 2 ist die gleiche Flä- 
che abgebildet, nachdem nochmals 
sorgfältig poliert worden war. 
Bringt man auf die so vorbereitete 
Probe einige Tropfen einer Fe,O,- 
Suspension, dann wird auch ohne 
Atzen das Gefüge. sichtbar, da sich 
über dem magnetischen Martensit 
eine stärkere Konzentration der Sus- 
pension einstellt als über dem un- 


Abb. 1. Tetragonaler und angelassener 
Martensit inmitten von Restaustenit auf 
Stahl; mit Salpetersäure geätzt. Die Härte- 
eindrücke dienen zur Markierung. V ~ 430 


magnetischen Restaustenit (Abb. 3). Wird ein Magnetfeld in der Größenord- 
nung von 200 Oe senkrecht zur Probenoberfläche angelegt, so zeigen sich auf 


1) Siehe z. B.: H. J. Williams, R. M. Bozorth u. W. Shockley, Physic. Rev. 75, 
155 (1949); R. M. Bozorth, J. Physique Radium 12, 308 (1951); L. F. Bates, Rev, 
Metallurg. 49, 185 (1952); B. Elschner, Ann. Physik 18, 290 (1953). 
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> La da 
Abb. 2. Dieselbe Stelle der Probenober- Abb. 3. Wie Abb. 2; mit einer Fe,0,-Sus- en! 
fläche wie in Abb. 1, jedoch nochmals pension bedeckt. (Durch das Zwischen- 


mechanisch poliert polieren hat sich die Form der Nadeln 
gegenüber der Abb. 1 etwas geändert) 


Abb. 4. Wie Abb. 3; ein Magnetfeld der Abb. 5. Wie Abb. 4; Magnetfeld um- 
Größenordnung von 200 Oe ist senkrecht gepolt = 
zur Oberfläche angelegt 
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W. Andrä: Magnetische Pulvermuster auf Martensit 


dem Martensit eigenartige Muster (Abb. 


33 
4). Beim Umpolen des Magnetfeldes 


vertauschen die hellen und dunklen Streifen ihre Lage (Abb. 5). 
Dieser Effekt, der in den Abb. 6—8 ebenfalls deutlich zu sehen ist, sowie 
die Form der Muster erinnern an die sogenannten Labyrinthmuster, die auf 


mechanisch poliertem  Siliciumeisen 
beobachtet werden können?). Ihre 
typische Form wird hierbei Spannun- 
gen zugeschrieben, die infolge des 
Schleifens und Polierens in der Ober- 
fläche entstehen. Die Labyrinthfiguren 
auf Silieiumeisen verschwinden nach 
einwandfreiem elektrolytischem Polie- 
ren, da dann die verspannte Schicht 
abgetragen worden ist. Der Charakter 
der Muster auf dem Martensit bleibt 
hingegen nach mechanischem und nach 
elektrolytischem Polieren der gleiche. 
Es kann daher angenommen werden, 
daß diese Muster durch magnetische 
Anisotropie im Inneren bestimmt sind, 
die im Martensit bei der Umwandlung 
entsteht und die sich natürlich nicht 
durch Abtragen einer Oberflächen- 
schicht beseitigen läßt. 

Wie sich aus dem Verhalten der 
Muster auf Martensit bei Veränderun- 


Abb. 6. Martensitnadeln auf sorgfältig 

mechanisch und elektrolytisch poliertem 

Stahl, mit einer Fe,0,-Suspension sicht- 
bar gemacht. V ~ 450 


gen des angelegten Magnetfeldes erkennen läßt, werden in diesem Fall 
keine Blochwände sichtbar gemacht. Es handelt sich vielmehr um eigen- 


Abb. 7. Wie Abb. 6; Magnetfeld senk- 
recht zur Oberfläche 


2) Siehe z. B.: L. W. McKeehan u. W. 


Abb. 8. Wie Abb. 7; Magnetfeld um- 
gepolt 


C. Elmore, Physic. Rev. 46, 226 (1934). 
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tümlich angeordnete Gebiete mit starker magnetischer Oberflächenladung. 
Da sehr ähnliche Muster auch bei Kobalt beobachtet werden konnten, liegt 
es nahe, dem Martensit ebenfalls eine einzige Richtung leichtester Magnetisier- 
barkeit und eine hohe Anisotropieenergie zuzuschreiben, wie N eel®) in an- 
derem Zusammenhang bereits ausgeführt hat. 

Weitere experimentelle Untersuchungen sowie Abschätzungen zur quan- 
titativen Deutung werden z. Z. durchgeführt. 


Jena, Institut für Magnetische Werkstoffe der Deutschen Akademie 
der Wissenschaften. 


Bei der Redaktion eingegangen am 17. April 1954. 


®) L. Neel, J. Physique Radium 12, 339 (1951). 
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Über die Phasenschwingungen beim Synchrotron 
mit starker Fokussierung 


(Eine Untersuchung mit Hilfe einer mechanischen Analogiemaschine) 


Von E. Bodenstedt 
Mit 17 Abbildungen 


Inhaltsübersicht 


Die Untersuchung zeigt, daß die mit der kritischen Energie verbundenen 
Erscheinungen keineswegs zu unüberwindlichen Schwierigkeiten führen. Es 
ist durchaus möglich, in der Beschleunigung von Teilchen durch die singuläre 
Stelle hindurchzukommen, ohne den Strahl zu verlieren, wenn man geeignete 
Eingriffe vornimmt. Es wird eine Anzahl soleher Maßnahmen näher disku- 
tiert, und gezeigt, daß durchaus ausreichende Toleranzen für die technische 
Verwirklichung vorhanden sind. 


Einführung 


Während beim normalen Synchrotron automatisch eine Stabilisierung der 
Phasenlage der Teilchen bezüglich der am Beschleunigungsspalt anliegenden 
Wechselspannung vorhanden ist, treten bekanntlich beim Synchrotron mit 
starker Fokussierung!) Schwierigkeiten auf, wenn man die Teilchen über die 
sogenannte kritische Energie hinweg beschleunigen muß. 

Im folgenden werden die eigentümlichen Verhältnisse in der Umgebung 
der kritischen Energie mit Hilfe einer mechanischen Analogiemaschine ein- 
gehend untersucht ; insbesondere wurden die Möglichkeiten diskutiert, die kri- 
tische Energie durch äußere Eingriffe zu überwinden. Es gelang, allgemein- 
gültige und von speziellen Synchrotronparametern unabhängige Aussagen 
darüber zu gewinnen, welche Amplitudenänderungen der Phasenschwingungen 
beim Durchgang durch die kritische Energie zu erwarten sind, und welche tech- 
nischen Toleranzen man bei den diskutierten Eingriffen einzuhalten hat. 


1. Über die der Untersuchung zugrunde liegende Differentialgleichung und 
ihren Gültigkeitsbereich 


D. Bohm und L. Foldy gaben eine Differentialgleichung für die Phasen- 
schwingungen beim normalen Synchrotron an?). Die folgende anschauliche und 


1) E. D. Courant, M. St. Livingstone and H. S. Snyder: The strong-focusing 
Synchrotron — A new high Energy-Accelerator, Physic. Rev. 88, 1190 (1952). 
as 2) D. Bohm and L. Foldy, The Theory of the Synchrotron, Physic. Rev. 70, 249 
46). 
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sich auf die wesentlichen Glieder®) beschränkende Herleitung läßt erkennen, 
wie die Gleichung auf den Fall der starken Fokussierung zu übertragen ist, 
und welche Vernachlässigungen hier den Gültigkeitsbereich einengen. 

Bekanntlich ist die Phasenschwingung ursächlich mit einer Energie- 
schwingung verknüpft. 

Die Phase ® eines Teilchens (Definition siehe Abb. 1) bestimmt die An- 
derung der Amplitude der Energieschwingung pro Umlauf A(AE) (AE = 
vn E—E, = Amplitude der Energieschwin- 

gung = Abweichung der Energie von 
der des synchronen Teilchens) : 


A (AE) = e- V- (sin — sin @,) *). (1) 


“=? Die Amplitude der Energieschwingung 
bestimmt aber riickwirkend die Ab- 
weichung der Winkelgeschwindigkeit 
Abb. 1. Zur Definition der Phase. V()= von der des synchronen Teilchens’) 
Beschleunigungsspannung; = Phase und damit die Anderungsgeschwindig- 
der Teilchen, definiert durch den Zeit- : x 
. keit der Phase: 
punkt, in dem sie den Beschleunigungs- 
spalt passieren; ®, = Gleichgewichtsphase 
(synchrone Phase), das Teilchen dieser 
Phase (synchrones Teilchen) erfährt pro 
Umlauf gerade den erforderlichen, durch : 
die zeitliche Zunahme des Führungsfeldee Diese von D. Bohm und L. Foldy 
und die Strahlungsverluste vorgegebenen verwandte Gleichung ist allerdings nur 
Energiebetrag eine Näherung, da es sich um eine 
Taylorentwicklung handelt, die nach 
dem ersten Glied abgebrochen wurde. Exakt gilt natürlich: 
do dE 
—=-—K. (2a) 


Die Vernachlässigung ist im Falle des normalen Synchrotrons ganz unwesent- 
lich, weil sich K nur langsam mit der Energie ändert, so daß die durch die 
Ableitungen von K bestimmten Koeffizienten der Glieder höherer Ordnung 
klein bleiben ; im Falle der starken Fokussierung macht sie sich jedoch in der 
unmittelbaren Umgebung der kritischen Energie (Vorzeichenwechsel von K) 
unter Umständen bemerkbar, wie am Ende dieses Abschnitts ausführlich 
gezeigt wird. 

Die Gl. (1) und (2) beschreiben den physikalischen Vorgang der Phasen- 
schwingung anschaulich und vollständig. Durch Elimination von AE gewinnt 
man unmittelbar die Phasenschwingungsgleichung von D. Bohm und 
L. Foldy: 


d/’E, 
= (a K' 6) + 5, (sin ®—sin®,) = 0. (3) 
3) Die Radialabhängigkeit der Ausstrahlung der auf einem Kreise umlaufenden Teil- 
chen führt zu einem von D. Bohm und L. Foldy angegebenen unwesentlichen Dämp- 
fungsglied, das hier vernachlässigt wurde. 
*) Das Teilchen der Phase ® gewinnt nämlich pro Umlauf den Energiebetrag 
e+ V-sin®, das synchrone Teilchen dagegen den Energiebetrag e- V - sin ®,. 
5) Der Momentanwert der Energie E bestimmt nämlich die Winkelgeschwindigkeit 


eindeutig; denn w = = , E legt unmittelbar die Bahngeschwindigkeit v fest und über 
das Führungsfeld auch den Radius r. 


B. Bo 
| wo & 
| ist, d 
fernt 
weic! 
(es is 
danr 
Eneı 
Insb 
dure 
(es i 
| d.h. 
1 
und 
Fok 
schv 
nur 
sein 
bek 
In 
eine 
Die 
füh 
beil 
nur 
Ap 
Neg wer 


E. Bodenstedt: Über die Ph i gen beim Synchrotron mit starker Fokussierung 37 


Allgemein lautet der Ausdruck für K: 
K=-1+@-1-5, (4) 


wo x eine, im allgemeinen von der Energie unabhängige Apparatekonstante 
ist, die angibt, um welchen Betrag sich ein Teilchen radial vom Sollkreis ent- 
fernt, wenn sein Impuls von dem des synehronen Teilchens p, um Ap ab- 
weicht : 
r ps 

Wenn x > 1 ist, dann behält X während der ganzen Beschleunigung sein 
Vorzeichen. Insbesondere ist beim normalen Synchrotron immer « > 1 
(es ist nämlich hier x = 1/1—n mit 0 < n< 1)’). Wenn dagegen x < 1 ist, 
dann wechselt K sein Vorzeichen, wie unmittelbar aus Gl. (4) folgt, bei der 
Energie : 


(4a) 


Exrit = m. ‚ wo E, = Ruhenergie. (5) 
& 


Insbesondere ist beim A.G.-Synchrotron (Synchrotron mit Fokussierung 
durch Sektorfelder mit alternierendem Feldgradienten)!) immer x < 1 
(es ist nämlich hier x = 4,85/n in asymptotischer Näherung für große n, 
d.h. etwa für n > 100)?). 

Es sei in dieser Arbeit ‚Starke Fokussierung‘ definiert durch « < 1 
und „Schwache Fokussierung‘ durch x > 1; physikalisch bedeutet „Starke 
Fokussierung‘ in diesem Sinne entsprechend Gl. (4a), daß große Impuls- 
schwankungen, was gleichbedeutend ist mit großen Energieschwankungen, mit 
nur kleinen Radialbewegungen verbunden sind. 

Da K in der Phasenschwingungsgleichung im Nenner auftritt, führt 
sein Vorzeichenwechsel zu einer unangenehmen Singularität. Der Term: 

E, 
K 
bekommt an der Stelle der kritischen Energie einen mathematischen Pol. 
In Abb. 2 ist die Funktion @(t) graphisch dargestellt unter der Voraussetzung 
einer linearen Zunahme der Energie des synchronen Teilchens mit der Zeit: 


E, = Exit (1 + k-t). (6) 


Alt) = 


Die Zeittransformation : 

t=k-t 
führt zu einem invarianten Zeitmaßstab, der in der folgenden Untersuchung 
beibehalten wurde. Beachtenswert ist das Minimum an der Stelle t = V 3—1. 


Es wurde oben erwähnt, daß die Phasenschwingungsgleichung von Bohm und Foldy 
nur eine Näherung ist. Sie ist im Falle des normalen Synchrotrons besonders gut und 

- folgt nämlich durch Differentiation: =n 
a) ‚AR _ 4a und unter Ver- 

p E 
{ c?) AE 
wendung der relativistischen Gleichung AE = v, - Ap: + («—1) 
7) n = Exponent der radialen Feldänderung: Bw 
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beschreibt auch im Falle des Protronensynchrotrons mit A.G.-Fokussierung selbst die 
unmittelbare Umgebung der kritischen Energie noch hinreichend genau, sie gibt dagegen 
im Falle des Elektronensynchrotrons mit A.G.-Fokussierung bei großen Amplituden und 
tiefer kritischer Energie die Umgebung der kritischen Energie nicht mehr ganz richtig 
; wieder. Exakter hat man an der Stelle von (2) die Gl. (2a) zu 
4 Olt, verwenden und gewinnt durch Integration: 


t Es+4 


E 
K(E) 


(2b) 


Es 


Die Elimination der Energie, die zur Phasenschwingungsgleichung 
; führt, macht jetzt Schwierigkeiten. Dagegen ist es ohne weitere 
4 , möglich, die Phase zu eliminieren, wenn man sich auf kleine 
= Phasenamplituden beschränkt, was in der unmittelbaren Um- 
gebung der kritischen Energie zulässig ist. Man kommt dann zu 
der Differentialgleichung der nit der Phasenschwingung gekop- 
pelten Energieschwingung. 


Unter Verwendung der 
Taylorentwicklung: 


sin = sin ®, + cos, (G—9) 
läßt sich nämlich Gl. (1) um 
formen zu: 


Man differenziert, um ® 
eliminieren zu können, einmal 
nach der Zeit: 


a? 
(AB) 


Die Zeitabhingigkeit vin o 
| und ®, wurde hier vernach- 
| Abb. 2. Die Funktion: ©(t) = E,/o} K. Zugrunde lässigt. Nun ist die Elimination 
i liegt ein linearer Anstieg der Energie des synchronen @ 1 d (2a) mög- 

Teilchens mit der Zeit: E,=E, „(I+kb lich, +) 


Es+4 
e-w-V-cos®, 
22 


AE + 


Unter Voraussetzung eines linearen Anstiegs der synchronen Energie mit der Zeit (6) 
kann man die Integration explizit durchführen und erhält, wenn man die synchrone 
Energie an Stelle der Zeit als unabhängige Variable einführt, unter Beachtung von Gl. (4) 
die folgende Differentialgleichung der Energieschwingung: 
2 


dy * 3 +3 

& = Energieamplitude AE/E, 

n = synchrone Energie E,/E,, 
Exrit = Kritische Energie E, 
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Interessant ist der Ort des Vorzeichenwechsels der Direktionskraft in cor &—n-Ebene. 
Diese Ortslinien mit dem Scharparameter ¢,,,, wurden errechnet und in Abb. 3 darge- 
stellt. Die Ortslinien haben bei kleinem £ einen Anstieg von ungefähr zwei. Das bedeutet 
insbesondere, daß ein Teilchen die kritische Stelle des Vorzeichenwechels der Direktions- 


kraftnicht dann erreicht,wenn#,=E, 


(definiert durch Gl. (5)), oder dann, 12 
wenn seine tatsächliche Momentanener- PEL OL Eee 


gieE = sondern dann, wenn 
gerade in der Mitte zwischen der tat- 
sächlichen und der synchronen Energie 
liegt. 


Das Diagramm ermöglicht folgen- h 
des Kriterium für die Gültigkeit 
der Phasenschwingungsgleichung von 
Bohm und Foldy in der Umgebung 
der kritischen Energie: 
-2r 


Notwendig und hinreichend ist 


| : ays Abb. 3. Die Ortslinien des Vorzeichenwechsels 
der Direktionskraft in der £ n-Ebene 


2. Eine mechanische Analogiemaschine zum Studium der Phasenschwingungs- 
gleichung 


Die Phasenschwingungen beim normalen Synchrotron und beim Syn- 
chrozyklotron wurden mehrfach in der Literatur durch Pendelmodelle veran- 
schaulicht?)®). In Weiterführung dieses Gedankens wurde eine mechanische 
Analogiemaschine entwickelt, die exakte Integrationen erlaubt und auch den 
Fall des Synchrotrons mit starker Fokussierung erfaßt. 


Es handelt sich also darum, ein mechanisches System zu finden, dessen 
Bewegungsgleichung mit der Phasenschwingungsgleichung identisch ist; 
denn man braucht dann nur den Bewegungsablauf zeitlich zu registrieren, um 
die Lösungen zu erhalten, Eine nähere Betrachtung der Phasenschwingungs- 
gleichung : 


(At) + (sin ® — sin ®,(t)) = 0 


zeigt, daß das folgende mechanische System diese Forderungen erfüllt: 


Man interpretiert ©(t) als das zeitlich veränderliche Trägheitsmoment 
eines Drehpendels (mit vertikaler Achse), ® als den Drehwinkel und damit ® 
als die Winkelgeschwindigkeit des Pendels. Das Produkt @(t) @P entspricht 
dann dem Drehimpuls. Die gesamte Gleichung wird zur mechanischen Be- 
wegungsgleichung des Pendels, wenn man das zweite Glied als ein Dreh- 
moment interpretiert, das auf das Pendel einwirkt. Die Gleichung schreibt 
also einmal ein Drehmoment vor, das dem Sinus des Ausschlagswinkels 


proportional ist, nämlich das Moment = 


schlagwinkel unabhängiges Moment 7 sin @,(t). Die schematische Skizze 


sin®, und ein zweites, vom Aus- 


8) M. Stanley Livingston, Particle Accelerators S. 291, Advances in Electronics 
Volume I, 1948, 
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auf Abb. 4 läßt erkennen, wie diese an das mechanische System gestellten 
Forderungen technisch verwirklicht wurden: 


Das Drehpendel hängt an einem Stahldraht von !/„ mm Durchmesser; 
der Durchmesser ist so klein gewählt, damit die Direktionskraft dieses Auf 
hängefadens gegenüber den übrigen Dreh- 
momenten vernachlässigt werden kann. Das 
Trägheitsmoment des Pendels wird durch 
zwei Bleimassen bestimmt, deren Abstand 
von der Drehachse variiert werden kann, 
Bei der gebauten Analogiemaschine wurde 
eine Anderung um einen Faktor 250 er 
reicht. Wie aus der schematischen Skizze 
ersichtlich, läßt sich die zeitliche Anderung 
des Trägheitsmoments durch einen nach 
unten in Richtung der Drehachse heraus- 
geführten Seilzug bedienen, auch während 
das Pendel schwingt, ohne dabei irgend- 
welche zusätzlichen störenden Drehmomente 
auf das Pendel auszuüben. Die durch 
die Differentialgleichung vorgeschriebenen 
Seilzug zur Steuerung des  Drehmomente wurden in folgender Weise 
hervorgerufen : 
Das Moment - sin ® wurde dadurch 

gewonnen, daß das Pendel mit einer kleinen, 
stromdurchflossenen Spule versehen wurde, die sich in einem äußeren hori- 
zontalen homogenen Magnetfeld befindet. Das übertragene Moment ist dann: 


D = Drehmoment 
D=M-H.sin® M = magnetisches Moment der Spule 
H = äußeres Magnetfeld. 


Man hat also die geforderte Abhängigkeit vom Ausschlagwinkel, und kann 
überdies eine eventuelle zeitliche Veränderung der Senderspannung V = V(t) 
durch Veränderung der Stromstärke in der Pendelspule und damit durch 
Veränderung von M verwirklichen. Das von ® und damit natürlich auch von 


V unabhängige?) Moment 5 sin ®,(t) wurde durch eine Spiralfeder ge- 


wonnen, die durch einige volle Umdrehungen vorgespannt wurde. Die Größe 
des übertragenen Drehmoments folgt zwar dem Hookesschen Gesetz, doch 
sind die tatsächlichen Pendelausschläge klein gegenüber diesem großen Vor- 
spannungswinkel, so daß das Moment praktisch vom Ausschlagwinkel ® 
unabhängig ist. Zur zeitlichen Steuerung des Moments läßt sich, wie aus der 
Skizze ersichtlich, die Spiralfeder an ihrem oberen Ende rasch auf- und ab- 
wickeln. Eine Schwierigkeit, die bei der Steuerung der zeitlich veränderlichen 
Größen auftritt, liegt in der Polstelle von ©(t). Einmal ist es unmöglich, das 


®) Das Produkt e- V -sin®, bedeutet die Energieaufnahme des synchronen Teil- 
chens pro Umlauf, die nur durch die zeitliche Zunahme des magnetischen Führungsfeldes 
und durch den Strahlungsverlust bestimmt ist, aber in keiner Weise durch die Sender- 
spannung. 
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Trägheitsmoment wirklich unendlich werden zu lassen. Statt dessen folgt 
man mit dem Trägheitsmoment der Funktion ©(t) bis zu einem sehr großen 
Wert unmittelbar vor der Polstelle; man läßt dann das Trägheitsmoment 
während des Durchgangs durch die Polstelle auf diesem Wert stehen, um es 
dahinter wieder mit ©(t) abnehmen zu lassen. Dadurch entsteht kein wesent- 
licher Fehler; denn nach dem Drehimpulssatz hat das in der Umgebung der 
Polstelle große Trägheitsmoment eine winzige Winkelgeschwindigkeit zur 
Folge, an der Polstelle selbst müßte momentan ® völlig verschwinden, während 
beim Modell die Winkelgeschwindigkeit auf einem äußerst kleinen Wert bleibt. 
Dies ist der einzige Unterschied. Hinter der Polstelle bewegt sich dann das 
Pendel mit völlig richtigem Drehimpuls 

weiter. 


Weiterhin läßt sich natürlich der 
Vorzeichenwechsel von @(t) nicht am 
Trägheitsmoment des Pendels vorneh- 
men. Statt dessen wechselt man das 
Vorzeichen der Drehmomente. 


Die Steuerung der zeitlich veränder- 
lichen Größen geschieht mit Hilfe von 
Programmscheiben. Sie steuern auch 
den Zeitpunkt des Vorzeichenwechsels 
der Drehmomente sowie den Startzeit- 
punkt des Pendels. Das Starten des 
Pendels mit kontrollierten Anfangs- 
bedingungen geschieht mit Hilfe einer 
magnetischen Zange, die das Pendel auf 
einem bekannten,eingestellten Anfangs- 
ausschlagwinkel festhält, bis sie durch 
die Programmscheibe geöffnet wird. 
Die Anfangswinkelgeschwindigkeit ist 
also immer null, 


Die Registrierung der Pendelbewe- 
gungen geschieht mit einem Lichtzeiger 
über einen Spiegel auf einen mit kon- 
stanter Geschwindigkeit bewegten Film. 
Zur Übersetzung der großen Pendel- 
ausschläge auf kleine Ausschläge des Spiegels wurde der Spiegel nicht un- 
mittelbar am Pendel, sondern am Aufhängefaden befestigt, dicht unterhalb 
des oberen Endes. 


Die Gesamtanordnung ist aus Abb. 5 ersichtlich. Pendel, Magnet und 
Registriervorrichtung sind auf einem Drehtisch montiert, dessen Drehachse 
mit der des Pendels zusammenfällt. Die Bedeutung dieses Drehtisches wird 
im letzten Abschnitt erläutert. Ein paar technische Daten mögen die Größen- 
verhältnisse wiedergeben : Die Feldstärke im Luftspalt betrug etwa 3000 Gauß. 
Die Amperewindungszahl der Pendelspule lag bei 20. Der volle Bereich der 
Änderung des Trägheitsmomentes entsprach einer Anderung der Schwingungs- 
zeit von etwa 0,6 sec bis auf 10 sec. 


Ann. Physik. 6. Folge, Bd. 15 4 


Abb. 5. Photo der Gesamtanordnung 
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3. Allgemeiner Überblick über die Phasenschwingungen, die damit gekoppelten 
Energieschwingungen und die besonderen Verhältnisse in der Umgebung der 
kritischen Energie 


In genügender Entfernung oberhalb der kritischen Energie werden für 
kleine Amplituden die Lösungen der Phasenschwingungsgleichung durch die 
von D. Bohm und L. Foldy angegebene?) adiabatische Näherung beschrieben: 


wo J und y Integrationskonstanten sind. 
In entsprechender Weise gilt in genügender Entfernung vor der kritischen 


Energie: 

20;-K t . {/e K cos ®,\} \. 
denn, wie aus der Phasenschwingungsgleichung unmittelbar ersichtlich, ist 
wegen des anderen Vorzeichens von © hier nicht mehr ®,, sondern 2 — ®, 
die stabile Gleichgewichtslage. Die schematische Skizze, Abb. 7a, zeigt die 
beiden Bereiche möglicher Phasenschwingungen vor und hinter der kritischen 
Energie. Diese Näherungen beschreiben die Bewegungen nur für hinreichend 
kleine Amplituden, da sie auf einer sich auf das lineare Glied beschränkenden 


Taylor-Entwicklung des Terms :_ (sin ® — sin ®,) der Phasenschwingungs- 


gleichung beruhen. 
Die Untersuchung mit der Analogiemaschine zeigte zunächst einmal, wie 
die Schwingungen bei großen Amplituden aussehen. In Abb. 6 ist eine Auf- 


Abb. 6. Photogramm einer Phasenschwingung großer Amplitude 


nahme einer Schwingung großer Amplitude wiedergegeben. Das Trägheits- 
moment des Pendels wurde bei dieserAufnahme konstant gelassen. Erwartungs- 
gemäß wird die Bewegung unter dem Einfluß des konstanten Drehmomentes 
unsymmetrisch. Außerdem wird die Schwingungsdauer jetzt amplituden- 
abhängig. Mit zunehmender Amplitude nimmt die Schwingungsdauer stetig 
zu. Die maximal mögliche Amplitude reicht bis zur labilen Gleichgewichts- 
lage, auf der das Teilchen theoretisch beliebig lange bleiben kann, was einer 
unendlichen Schwingungsdauer entspricht. 

Weiterhin zeigte die Untersuchung bei kleinen Amplituden, daß der 
Bereich, in dem merkliche Abweichungen von der adiabatischen Näherung 
auftreten, sich auf die unmittelbare Umgebung der kritischen Energie be- 
schränkt, und zwar auf eine Zone von etwa zehn vollen Phasenschwingungen. 
In der adiabatischen Näherung gehen die Amplituden beim Annähern an 
die kritische Energie gegen null. In der exakten Lösung nehmen dieAmplituden 
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zwar auch beim Annähern an die kritische Energie stetig ab, doch bleiben 
sie an der kritischen Stelle selbst endlich. In Abb. 7 ist eine Aufnahme der 
Phasenschwingungen in der unmittelbaren Umgebung der kritischen Energie 
wiedergegeben. Das Teilchen schwingt zunächst um die Gleichgewichtslage 
z—®, herum. Frequenz und Amplitude nehmen ab. An der Stelle der 
kritischen Energie wird x — ®, zur labilen und ®, zur stabilen EP a 
lage. Das Teilchen befindet sich gerade etwas unterhalb von 2 — . (Die 
Stelle der kritischen Energie ist daran zu erkennen, daß dort immer = 6) 
Es führt jetzt eine weitere Schwingung um die neue Gleichgewichtslage ®, 


2 
— stabiler Bereich vor Exrjy — \— stabiler Bereich hinter Eqn, — 


Abb. 7. Photogramm, Ver- Abb. 7a. ~~ beiden Bereiche stabiler 
lust eines Teilchens hinter Schwingungen 
der kritischen Energie 


herum aus mit beinahe der maximal möglichen Amplitude. Hinter der kriti- 
schen Energie nehmen aber die Phasenschwingungsamplituden wieder zu, 
so daß das Teilchen über die labile Gleichgewichtslage x—®, hinausschwingt !P). 
Es verläßt jetzt den Bereich stabiler Schwingungen, seine Phase nimmt 
weiter beschleunigt zu (im Modell unter der Wirkung des konstanten Drek- 
momentes), und seine Energie sinkt immer tiefer unter den Wert der Energie 
des synchronen Teilchens, bis es auf die Innenwand der Beschleunigungs- 
kammer auftrifft und verloren geht. Um die Umgebung der kritischen Energie, 
insbesondere das Verhalten unmittelbar oberhalb der kritischen Energie ein- 
gehender zu studieren, wurde ein ,,Verlust der Teilchen‘‘ dadurch vermieden, 
daß kein Vorzeichenwechsel an den Momenten vorgenommen wurde. Die 
auf diese Weise erhaltenen Teilchenbahnen haben keineswegs nur theoreti- 
sches Interesse. Sie stellen nämlich genau die Bewegung dar, die beim Syn- 
chrotron auftreten würden, wenn man genau im Zeitpunkt des Erreichens der 
kritischen Energie momentan die Phase des Senders um den Betrag 2 — 2 ®, 
verschieben würde. 

Dieser Idealfall läßt sich zwar bei einem wirklichen Synchrotron niemals 
exakt verwirklichen, doch sind für eine angenäherte, mit technischen Mitteln 
mögliche Verwirklichung ausreichende Toleranzen da, wie im letzten Ab- 
sehnitt ausführlich beschrieben wird. Die Abb. 8 zeigt die Aufnahme einer 
Phasenschwingung, bei der der Durchgang durch Exrit in der beschriebenen 
Weise mit der Analogiemaschine vorgenommen wurde. 

Es möge an dieser Stelle eine Überlegung wiedergegeben werden, die die 
Verallgemeinerung spezieller Resultate auf beliebige Parameter betrifft. 


1) Vergleiche dazu die schematische Skizze auf Abb. 7a. 
4* 
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Die Phasenschwingungsgleichung enthält drei wesentliche Synchrotron- 
parameter: der erste wesentliche Parameter gibt an, wie rasch die synchrone 
Energie mit der Zeit zunimmt. Durch die oben erwähnte Zeittransformation 
t =k-t wurde er aus allen Überlegungen und Resultaten eliminiert. Der 

y zweite wesentliche Parameter 

ist die synchrone Phase 9, ; sie 

kann zeitabhangig sein. In dem 
kurzen Zeitintervall der Um- 
gebung der kritischen Energie, 
das im folgenden interessiert, 
möge sie jedoch als konstant 
angesehen werden. Sie be- 
stimmt dıe Gleichgewichtslage, 
um die herum die Schwingun- 
gen auftreten. Durchdie Trans- 
formation: ¥ = ®—@, hinter 
der kritischen Energie, bzw. 
WY = ®—(x—®,) vor der kri- 
tischen Energie läßt sich da- 
mit auch @, eliminieren, Tat- 
sichlich hat man im allgemei- 
nen in der Umgebung der 
kritischen Energie wegen der 
vorangehenden Amplituden- 
abnahme mit kleinen Ampli- 
tuden zu rechnen, so daß 
man schon aus diesem Grunde 
zweckmäßigerweise die Unter- 
suchung auf kleine Amplitu- 
den beschränkt. Der dritte 
wesentliche Parameter ist die 
Stärke der Phasenfokussie- 
rung bzw. die Direktionskraft 
in der Phasenschwingungsglei- 
ehung. Sie ist im wesentlichen 
durch die Senderspannung be- 

Abb. 8. Durchgang durch die kritische Energie mit Stimmt (und, wenn man nicht 
Hilfe des idealen Phasensprungs. Darstellung von: mit der Grundschwingung ar- 
Y(t), @(t) und AE beitet, sondern mit einer Ober- 

welle, durch deren Ordnung); 

einen kleinen Einfluß auf die Direktionskraft hat außerden ®,. Die Stärke der 
Phasenfokussierung bestimmt natürlich die Phasenschwingungsfrequenz und 
damit die zeitliche Dauer der interessierenden etwa zehn vollen Phasen- 
schwingungen in der Umgebung von Exrit, die durch die adiabatische Nä- 
herung nicht beschrieben werden. Diese Zeitspanne sei: „<r<+tYn 
(7, repräsentiert somit den Parameter der Stärke der Phasenfokussierung). 
Dieser dritte Parameter hat im allgemeinen einen grundsätzlichen Einfluß 
auf die Gestalt der Lösungen, da er das durchlaufene Programm der Funktion 
Or) bestimmt. Allerdings ändert sich das Programm, wenn man mit 7, 
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heruntergeht, unterhalb von etwa t, = 0,05 praktisch verhältnisgleich, so 
daß unterhalb von r, = 0,05 die Lösungen auch gegenüber t, invariant werden. 
Aus diesem Grunde wurden die folgenden Messungen ausschließlich bei 
%< 0,05 ausgeführt. Die Bevorziehung dieses Gebietes ist überdies deswegen 
gerechtfertigt, weil man bei é 

großen Protonensynchro- 
trons mit starker Fokussie- 
rung wohl auch wirklich 
in diesem Gebiet arbeiten 
wird. 

Die Lésungen der Pha- 
senschwingungsgleichung 
unterscheiden sich schlieB- 
lich noch durch ‚zwei Inte- 
grationskonstanten, alsoet- 
wa Amplitude und Phase. 
Eine Variation der Ampli- 
tuden führt zu mathema- 
tisch ähnlichen Lösungen, 
solange die Amplituden 
klein sind, da, wie unmit- 
telbar aus der Gestalt der 

Phasenschwingungsglei- 
chung folgt, dann mit Y(t) - 
auch c-Y(r) eine Lösung 
ist. Die Phasenlage der 
Phasenschwingungdagegen 
bleibt ein wesentlicher 
Parameter, der sorgfältige 
Berücksichtigung verlangt. 
Im Sinne dieser Überlegun- 
gen ist die in Abb. 8 Apb.9 Enveloppen für ® und AE für kleine Ampli- 
dargestellte Phasenschwin- tuden der Phasenschwingung 
gung eine durchaus allge- 
meingültige Lösung der Phasenschwingungsgleichung in dem Bereich, in dem 
die adiabatische Näherung versagt. 

Es muß an dieser Stelle eine Untersuchung von K. Johnson!!) erwähnt 
werden, dem es gelang, von derselben Phasenschwingungsgleichung ausgehend 
für die unmittelbare Umgebung der kritischen Energie eine analytische 
Näherungslösung zu gewinnen. Auf Abb. 8 ist unter der Phasenschwingung 
ihre Ableitung aufgetragen und darunter die mit der Phasenschwingung ge- 
koppelte Energieschwingung. Man gewinnt die Energieschwingung aus der 
Phasenschwingung unter Verwendung der Gleichung: 


11) K. Johnson, Phase Oscillations and Transition Energy Problems, Lectures on 
the Theory and Design of an Alternating-Gradient Proton Synchrotron (Conference of 
the Alternating-Gradient Proton Synchrotron at Geneva, Oktober 26-27-28, 1953. 
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—kt=T 


Abb. 10. Die 4 Typen des Durchgangs durch die kritische Energie beim idealen 
Phasensprung des Senders 
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oder, da: Aw = ® 
AE = Er 
= - O(t)-O-w,. 


Man muß also die Phaserschwingung nach der Zeit differenzieren und im 
wesentlichen mit dem Trägheitsmoment ©(t) multiplizieren. Für O und AE 
gibt es Enveloppen ; sie wurden mit der Analogiemaschine aufgenommen und 
in Abb. 9 dargestellt. Interessant ist der Pol der Enveloppe der Energie- 
schwingung an der Stelle der kritischen Energie. Es ist wesentlich für das Bild, 
das man sich von den Vorgängen an der kritischen Stelle zu machen hat, daß 
die Amplitude der Energieschwingung beim Annähern an die kritische Energie 
zunimmt. Sie kann zwar tatsächlich nicht unendlich werden, doch schwingt 
die Energie zwischen den Enveloppen an der Polstelle um so höher hinauf, je 
größer die Phasenschwin- 


gungsfrequenz, d.h. aber je 
stärker die Phasenfokussie- 
rung ist. Dies ist insofern 
von Bedeutung, als mit der 
Energieschwingung natürlich 
eine Radialbewegung gekop- 
pelt ist, deren Amplitude 
nicht größer werden darf als 
der Radius der Kammer. 0 1/2 1 Jor Zr 


Interessant ist folgendes App. 11. Die Maximalamplitude ¥,,, an der Stelle 


weitere Resultat für den t = + 7, als Funktion der Phasenlage der Phasen- 
Durchgang ‚durch die kriti- schwingung p an der kritischen Energie. Prax an 
sche Energie mit Hilfe des der Stelle r = — Tu beträgt 22° und es liegen volle 
„Idealen Phasensprunges‘‘. 10 Phasenschwingungen dazwischen (gültig für 


Die Ordinaten der Enveloppe 7» < 0,1). Zur Definition von 9: 


an die Phasenschwingung 

sind in symmetrischen Punk- SN aaa 

ten auf der Zeitskala vor und N >t 
hinter der kritischen Energie ee i i» 

selbst in diesem idealen Fall ” Zu 


keinesfalls gleich, sondern 

die Amplitude ändert sich um einen Faktor, der abhängig ist von der 
Phasenlage der Phasenschwingung, in der die Teilchen die kritische Stelle 
erreichen. Abb. 10 zeigt die vier typischen Fälle. Die beiden ersten Fälle 
sind die Grenzfälle des symmetrischen und antisymmetrischen Durchgangs, 
bei denen sich die Amplituden nicht ändern. Die ersichtliche Verschiebung 
der Phasenlage zwischen Phasenschwingung und Energieschwingung hat ihre 
Ursache natürlich in der Differentiation nach der Zeit. Die beiden anderen 
Fälle sind Beispiele für Zunahme und Abnahme der Amplitude. Der Verlauf 
der Schwingungen selbst erlaubt es, anschaulich zu verstehen, wie eine 
Anderung der Amplituden zustandekommt. Die Abb. 11 zeigt den Faktor 
der Amplitudenänderung zwischen symmetrischen Punkten als Funktion der 
Phasenlage der Phasenschwingung an der kritischen Stelle. 
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4. Über die technischen Möglichkeiten, die kritische Energie zu überwinden 
a) Ein Phasensprung des Senders 


Im Zeitpunkt des Erreichens der kritischen Energie läßt man die Phase 
des Senders um den Betrag 7 —2®, springen. Diese Möglichkeit wurde 
bereits in der Arbeit von E. D. Courant, M. S. Livingston und H.$, 
Snyder erwähnt!). 

Ein solcher Phasensprung kann technisch jedoch nur in endlicher Zeit 
ausgeführt werden. Außerdem kann er nur für wenige Teilchen angenähert 
im richtigen Zeitpunkt erfolgen ; die tatsächlich im Synchrotron schwingenden 
Teilchen erreichen nämlich die kritische Stelle des Vorzeichenwechsels der 
Direktionskraft zu verschiedenen Zeitpunkten aus zwei Gründen: 

Der erste Grund ist in den Überlegungen am Ende des ersten Abschnitts 
enthalten. Er liegt darin, daß die synchrone Energie nicht allein den Zeit- 
punkt des Vorzeichenwechsels der Direktionskraft bestimmt, sondern daß 
auch die momentane von Teilchen zu Teilchen verschiedene Energie- 
schwingungsamplitude einen Einfluß hat, wie in Abb. 3 dargestellt wurde, 

Ein zweiter Grund, auf den K. Johnson") aufmerksam machte, liegt 
insbesondere beim A. G. Synchrotron in Nichtlinearitäten der radialen Feld- 
änderung. Ist nämlich » radialabhängig, so wird damit auch die von n ab- 
hängige kritische Energie zu einer Funktion des Radius, und da die Teilchen 
sich wegen der mit der Phasenschwingung gekoppelten Radialschwingung auf 
verschiedenen Mometanradien bewegen, unterscheiden sie sich damit auch im 
Wert von Exrit- Das Ziel der folgenden Messungen war, allgemeingültige 
Aussagen über die Toleranzen für die Dauer des Phasensprungs sowie für 
den zeitlichen Abstand zwischen Phasensprung des Senders und Vorzeichen- 
wechsel der Direktionskraft beim Teilchen zu gewinnen. 

Ein Phasensprung des Senders in endlicher Zeit läßt sich bei der Analogiema- 
schine in besonders einfacher Weise dadurch verwirklichen, daß man den Dreh- 
tisch in dieser Zeit um den Winkel der gewünschten Phasen verschiebung dreht. 

Ein Phasensprung des Senders um den Wert Ag in der Zeit von t, bis t, ist nämlich 
äquivalent einem kurzzeitigen Frequenzfehler Aw mit: 

tz 

Ag = J Aw - dt. 

tı 
Der Frequenzfehler Aw läßt sich gerade dadurch verwirklichen, daß man den Drehtisch 
mit der Winkelgeschwindigkeit Aw dreht: Die Phasenschwingungsgleichung beschreibt 
die Schwingung zwar für beliebige Zeitabhängigkeit von w,. Aber kleine Schwankungen 
um den richtigen Frequenzgang herum führen schon zu starken Schwankungen von sin ®,, 
die unbequem zu berechnen und noch schwieriger zu steuern wären. Es ist sehr viel 
richtiger, die für den idealen Frequenzgang berechneten Funktionen ©(t) und sin ®,(t) 
zu verwenden und statt dessen den Einfluß von Frequenzfehlern durch eine Taylor- 
entwicklung zu gewinnen. Diese sei im folgenden durchgeführt. 

Es sei also: 

o,= + Aw 
und damit: 


+ Aw (t) = momentaner Frequenzfehler. 
dAw _ dos 


Es sei: Aw < @,9, jedoch keineswegs: In der Phasenschwingungs- 


©, = tatsächlicher Frequenzgang 
sg = idealer Frequenzgang 


gleichung (3) sind die Glieder: O(t) = are » sowie — - sin ®,(t) zu entwickeln. Die 
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Entwicklung von ©(t) liefert: 
u 3 2 2\ Aw Es 
tt 
Alle Glieder klein von erster Ordnung mögen vernachlässigt werden, da sie nichts aus- 
machen im Vergleich zu dem Hauptglied, das bei der Entwicklung von sin ®,(t) auftritt: 
L, 
an sin ®, = + . 
Diese Gleichung ist nichts anderes als die Energiebilanz. (Z, = Strahlungsverlust/Um- 
lauf.) Mit: E,= E,, + AE ergibt sich bei Vernachlässigung der Glieder klein von erster 
Ordnung: 


(Es unterscheidet sich nämlich L, nur um ein Glied klein von erster Ordnung von Li.) 
Es ist AZ keineswegs klein von erster Ordnung gegen E,,, da es nicht durch Aw, sondern 


durch ze bestimmt wird. 
Es tritt also das Zusatzglied auf: 


ey. e 
Dies Zusatzglied ist jetzt mit Hilfe des Frequenzfehlers auszudrücken, was unter Ver- 
wendung der Gl. (2): 
Aw, AE 


so E 30 


gelingt. Nach Differentiation folgt nämlich: 


Damit lautet die neue Differentialgleichung: 
dr Es ER; 
Fas (Ö+ Aw) | + 5— (sin ® — sind) = 0. 


Sie unterscheidet sich von der alten Gleichung nur dadurch, daß zu ® der Frequenz- 
fehler additiv hinzutritt. Dies additive Zusatzglied läßt sich im Modell nun aber gerade 
dadurch verwirklichen, daß man den 
ganzen Mechanismus auf einen Dreh- 
tisch setzt und mit der Winkelgeschwin- 
digkeit Aw rotieren läßt. Die neue 
Bewegungsgleichung des Pendels für 
das Koordinatensystem des mitbeweg- 
ten Beobachters (mitbewegte Kamera) 
ist dann nämlich mit dieser neuen Glei- 
chung identisch, wie leicht einzusehen 


“te Abb. 12 sind zwei Photo- 


gramme von Teilchenbahnen wie- Abb, 12. Photogramme der Phasenschwingung 
dergegeben, bei denen der Phasen- beim Durchgang durch E, . mit Hilfe eines 
sprung des Senders in der be- Phasensprungs von endlicher Dauer 
schriebenen Weise vorgenommen 

wurde. Das Teilchen auf der ersten Aufnahme bleibt fokussiert. Bei der 
zweiten Aufnahme kam der Phasensprung nicht genau im richtigen Zeit- 
punkt. Das Teilchen hat nach dem Passieren von Exrit eine stark vergrößerte 


Amplitude und geht nach wenigen Schwingungen verloren. 
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In einer ersten Meßreihe wurde untersucht, wie sich die Amplituden 
zwischen symmetrischen Punkten (auf der Zeitskala) vor und hinter der 
kritischen Energie ändern als Funktion der Dauer des Phasensprungs. Der 
Phasensprung kam bei dieser Meßreihe also genau im richtigen Zeitpunkt. 
Da die Anderung der Amplituden von der Phasenlage der Phasenschwingung 
abhängig ist, mußte diese natürlich variiert werden. Abb. 13a zeigt das Resul- 
tat für eine Anfangsamplitude von 20° an der Stelle t = — r,. (tT) wie oben 
definiert dadurch, daß zwischen —r, und +7, zehn volle Phasenschwingungen 


Wo, 


] 


Abb. 13a. %, „, an der Stelle t= + T, 


als Funktion der Phasenlage der Phasen- 
schwingung an der Stelle der kritischen 
Energie. Scharparameter ist die Dauer 


Abb. 13b. Die Bereiche der Maximal- 
amplitude „, an der Stelle t= + % 


als Funktion der Dauer des Phasen- 
sprungs des Senders. //// Fan (t= 


ax 
(v7 = — = 0° 


des Phasensprungs des Senders: — = 20°,-+- 
m 
Kurve 


ax 


Dauer des Phasensprungs 

ind. Einh. ty 0 | 0,2) 0,44 

YP max an der Stelle r = — r, betrug 20° 

und zwischen —r, und +7, lagen 10 

volle Phasenschwingungen. Giiltig fiir 
alle t, < 0,1 


liegen). Das Resultat dieser Untersuchung läßt sich nicht unmittelbar für 
beliebige Amplituden verallgemeinern. In Abb. 13b sind die Bereiche der 
hinter der kritischen Energie an der Stelle r = + r, auftretenden Amplituden 
noch einmal für die Anfangsamplitude von 20° sowie für die Anfangsamplitude 
von 0° als Funktion der Dauer des Phasensprungs dargestellt. Insbesondere 
entnimmt man diesem Diagramm, daß der Faktor der Amplitudenänderung 
bei gleicher absoluter Dauer (in sec) des Phasensprungs um so größer wird, je 
stärker man die Phasenfokussierung macht; starke Phasenfokussierung 
heißt, wie in Abschnitt 3 gezeigt, kleines t,, und das bedeutet in Abb. 13b 
einmal eine große Abszisse, zum andern aber eine kleine Anfangsamplitude. 
Durch beides wird der Faktor der Amplitudenänderung groß. 


In einer zweiten Meßreihe wurde der Faktor der Amplitudenänderung 
untersucht, die auftritt, wenn der Phasensprung zwar momentan erfolgt; 
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doch nicht im richtigen Zeitpunkt. In diesem Fall spielt die Phasenlage der 
Phasenschwingung eine besonders starke Rolle. Die Resultate sind in den 
Abb. 14a und 14b gestrichelt eingetragen. Insbesondere erlauben die Kurven 
der Abb. 14a die Belegungswahrscheinlichkeit des Bereichs der Amplituden- 
faktoren, der in Abb. 14b dargestellt ist, zu berechnen. Die Größenordnung 
des Effektes läßt sich dadurch charakterisieren, daß die Amplituden beim 
Durchgang durch die kritische Energie im Mittel etwa um einen Faktor zehn 
zunehmen, wenn zwischen Phasensprung des Senders und dem Erreichen der 
kritischen Stelle durch das Teilchen eine halbe Phasenschwingung vergeht. 


Ya 


Abb. 14a. Die Abhängigkeit des Faktors der Amplitudenänderung von der Phasen- 
lage der Phasenschwingung 
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Abb. 14b. Die Bereiche des Faktors der Amplitudenänderung zwischen symmetrischen 
Punkten vor und hinter der kritischen Energie als Funktion des Zeitpunkts, an dem die 
krit. Stelle erreicht wird. \\\\ bei einem Phasensprung des Senders; //// bei 11 Phasen- 
springen des Senders (bei 10 vollen Phasenschwingungen zwischen —t, und +7, und 
bei kleinen Amplituden der Phasenschwingung). Das Diagramm ist gültig für 7, << 0,05 
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b) Eine größere ungerade Zahl von äquidistanten Phasensprüngen des Senders) 


Eine weitere Methode wurde in einer früheren Notiz erwähnt), Es ist 
die Anwendung des Prinzips der A. G.-Fokussierung (Sektorfelder mit alter. 
nierendem Feldgradienten) auf die Phasenschwingungen. Der Nachteil der 
Methode in der dort beschriebenen Form liegt darin, daß diePhasenfokussierung 
im allgemeinen hierbei sehr schwach werden dürfte, wenn man nicht zur 
Beschleunigung der Teilchen eine sehr hohe Spannung und eine Oberwell 
sehr hoher Ordnung verwendet. Der besondere Vorteil ist, daß alle Teilchen 
fokussiert bleiben unabhängig davon, wann sie die kritische Stelle erreichen. 
Man kann den erwähnten Nachteil durch folgende veränderte Anwendung de 
Prinzips vermeiden. Man arbeitet zunächst mit normaler Phrasenfokussierung. 
Kurz vor der kritischen Energie beginnend macht man nun in gleichen Zeit- 
abständen einige Phasensprünge des Senders um 2 — 2 ®, hin und her und 
erzielt auf diese Weise nach dem neuen Fokussierungsprinzip weiterhin 
Phasenfokussierung. Den Zeitabstand wählt man so, daß man etwa in der 
Mitte des Stabilitätsdiagramms bleibt. Dann ist die Stabilisierung nicht viel 
schlechter als vorher (nur etwa um einen Faktor zwei). Innerhalb dieses 
Zeitraums der Phasensprünge erreicht dann das Teilchen die kritische Energie 
und bleibt weiterhin fokussiert, fokussierender und defokussierender Zeit- 
abschnitt vertauschen lediglich ihre Rolle. Nach einer ungeraden Gesamtzahl 
von Phasensprüngen läßt man dann die Senderphase wieder konstant und 
hat wieder normale Fokussierung. 


Es wurde eine Meßreihe durchgeführt, bei der elf Phasensprünge des Senders 
verwandt wurden. Untersucht wurde wie vorher bei nur einem Phasen- 
sprung der Amplitudenänderungsfaktor als Funktion der Zeitdifferenz zwischen 
dem mittleren (d.h. dem 6.) Phasensprung und dem Erreichen der kritischen 
Stelle. Die Resultate wurden in den Abb. 14a und 14b ausgezogen eingetragen. 
Die Abb. 14b erlaubt einen aufschlußreichen Vergleich mit den Verhältnissen 
bei nur einem Phasensprung. Man erkennt, daß jetzt die Faktoren der Ampli- 
tudenänderung klein bleiben, solange das Teilchen die kritische Stelle nur 
irgendwann in dem gesamten Bereich der elf Phasenspriinge erreicht. Erst 
außerhalb der Zone der Phasensprünge steigt der Bereich der Amplituden- 
faktoren an, und zwar dort etwa genau so steil, wie er im Falle nur eines Phasen- 
sprungs in dessen Umgebung ansteigt. 


Wenn also der Streubereich für die Zeitpunkte, in denen die einzelnen 
Teilchen die kritische Stelle erreichen, groß ist, d. h. größer als die Dauer einer 
Phasenschwingung, dann ist die Methode eines Phasensprungs ungeeignet, 
dagegen die Methode einer größeren ungeraden Zahl von Phasensprüngen die 
ideale Lösung. 


¢) Modulation der Senderspannung 
Die schematische Skizze Abb. 7a legt folgende weitere Möglichkeit nahe: 
Die beiden Zonen stabiler Schwingungen vor und hinter der kritischen Energie 
überlappen sich praktisch nicht wegen der starken Dämpfung der Amplituden 


12) Die Anregung zu dieser Untersuchung verdanke ich Herrn Prof. D. Kerst, Illinois. 
13) E. Bodenstedt, Zur Überwindung der kritischen Energie beim Synchrotron mit 
starker Fokussierung, Z. Naturforschg. 8a, 502 (1953). 
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vorher und der Aufblähung dahinter. Es ist jedoch leicht möglich, die 
beiden Gleichgewichtslagen x — ®, und ®, an der kritischer. Stelle soweit zu 
nähern, daß nun doch eine $ 

Uberlappung der Zonen stabi- 
ler Schwingungen zu maximal 
50% eintritt. Die Steuerung 
yon ®, kann man dadurch vor- 
nehmen, daß man die Sender- 
spannung moduliert. Die tech- 
nische Durchführung dürfte 
gar keine Schwierigkeiten ma- 
chen, da diese Modulation nur 
ganz langsam vor sich zu gehen 
braucht. Schwierigkeiten dürf- mn 
ten allerdings in jedem Fall fiir stabiler Bereich vor Erit. — Bech hint 


die Anwendung dieser Methode 
bestehen, wenn man es an der Abb. 15. Schematische Darstellung der Uberlappung 
der beiden stabilen Bereiche durch Spannungs- 
kritischen Stelle mit besonder S modulation (diese ermöglicht eine Steuerung von ®, 
kleinen Phasenamplituden zu und damit auch von 2 — ®,) 
tun hat, wie es praktisch 
bei großen Protonenbeschleu- 
nigern der Fall sein dürfte. 
Abb. 15 veranschaulicht die 
Methode in einem schemati- 
schen Diagramm. Abb. 16 ist 
ein Photogramm, das mit der 
Analogiemaschine gewonnen 
wurde, und das die Überwin- 
dung der kritischen Energie 
mit der Methode der Span- Abb. 16. Photogramm des Durchgangs durch 
nungsmodulation wiedergibt, Ex... durch Spannungsmodulation 


t 


d) Die Toleranzen fiir die Steuerung der Senderfrequenz 


Eine mit der Überwindung der kritischen Energie verbundene ernste tech- 
nische Schwierigkeit liegt in den Toleranzen, die man bei der Steuerung der 
Senderfrequenz einzuhalten hat. Da die Teilchen beim Durchgang durch die 
kritische Energie die Lichtgeschwindigkeit noch nicht erreicht haben, nimmt 
ihre Bahngeschwindigkeit und damit ihre Winkelgeschwindigkeit weiterhin 
zu. In entsprechender Weise muß die Senderfrequenz zunehmen. Wegen der 
großen Dauer einer Phasenschwingung in der kritischen Zone ist die Bewegung 
in dieser Zone besonders empfindlich gegenüber Fehlern in der zeitlichen 
Steuerung der Senderfrequenz; das entspricht bei der Analogiemaschine also 
Bewegungen des Drehtisches. Ein zeitlich konstanter Frequenzfehler Aw 
verschiebt den Zeitpunkt, an dem die kritische Energie erreicht wird und den 
Sollkreisradius, hat aber auf die Phasenschwingung selbst sonst keinen weiteren 
Einfluß. Dagegen reagieren die Phasenschwingungen empfindlich auf zeitliche 
Änderungen des Frequenzfehlers; d.h. also auf Fehler in der Richtung der 
Frequenzzeitkurve. Ein solcher Fehler läßt sich mit der Analogiemaschine 
sehr einfach verwirklichen. Man braucht nur den Drehtisch konstant zu be- 
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schleunigen. Die Resultate einer Meßreihe sind in Abb. 17 wiedergegeben. 


Eine parameterfreie Darstellung war dadurch möglich, daß die zeitliche Anderung 


des Senderfrequenzfehlers 


2 
ante Zeiteinheit r, transformiert wurde durch den Faktor = . Man entnimmt 


von der Zeiteinheit Sekunde auf die invari- 


dem Diagramm zunächst einmal, daß die oberhalb der kritischen Energie auf- 
tretenden Amplituden Pax (tT = T)) nur wenig von der Anfangsamplitude 


Wnar(TTo) 


Abb. 17. Toleranzen für die Exaktheit des Frequenzhubs. Die Maximalamplitude an 

der Stelle r = r, als Funktion des Fehlers im Senderfrequenzhub (zwischen t = — y 

und t= +r, liegen 10 Phasenschwingungen, t)< 0,1). //// max (t = To) = 10°; 

—+—P ax 0°; w, = Kreisfrequenz des Senders; Aw, = Momentanwert 
der Abweichung der Kreisfrequenz des Senders vom Sollwert 


max (t = — r,) abhängen. Man ersieht weiter, daß es günstig ist, mit mög- 
lichst starker Phasenfokussierung zu arbeiten, da Pax (t = + r,) bei kon- 


stantem Fehler im Senderfrequenzhub = etwa mit 7% zunimmt, während 


die zulässigen Amplituden !4) nur mit Y% ansteigen, wie sich ohne besondere 
Mühe aus der adiabatischen Näherung gewinnen läßt. Durch Einsetzen 
spezieller Synchrotrondaten für 7, und k entnimmt man dem Diagramm 
unmittelbar die Größe des Effektes. 


Ich möchte an dieser Stelle Herrn Prof. Dr. W. Paul für die Förderung 
der Arbeit meinen herzlichen Dank aussprechen. Ich verdanke Herrn Werner 
Hardt und meinen Kollegen von der Synchrotrongruppe Dr. Ehren berg 
und Steinwedel zahlreiche wertvolle Diskussionen. Die Untersuchung 
wurde mit Mitteln der Deutschen Forschungsgemeinschaft durchgeführt. 


14) Sie müssen hinreichend klein sein, um bei der stetigen Amplitudenzunahme bis 
zu der Stelle des Minimums von @ nicht über den stabilen Bereich hinauszugehen. 


Bonn, Physikalisches Institut der Universität. 


Bei der Redaktion eingegangen am 28. April 1954. 
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Berichtigung zur Arbeit: 


„Einige Bemerkungen zur Gisolfschen Theorie der 
Elektronenschwankungserscheinungen von Halbleitern“ 


Ann. Physik 6 14, 87, 1954 
Von K. W. Böer 


Bei der Berücksichtigung statistisch verteilter Lebensdauern von Lei- 
tungselektronen wurde vom Verf. seinerzeit übersehen, daß die Elektronen- 
konzentration, welche die Zahl der sekundlichen Anregungsprozesse N und 
die Lebensdauer r der Leitungselektronen enthält, nicht vor das Integral in 
Gl. (4) gezogen werden darf. 


Unter Berücksichtigung dieser Tatsache und nach Verwendung von 
n= N r,/Lg ergibt sich an Stelle von (4) 


oo 

| sin? dr. (4*) 

Nach Ausführung der Integration erhält man statt (5) die mit den Ergeb- 

nissen von Bernamont und Van der Ziel übereinstimmende Gleichung 

4ebEr, Av 


R 1+ w 73° 


(5*) 


2 
„= 


In gleicher Weise sind auch im folgenden Paragraphen 22 die Integra- 
tionen (7) bzw. (7a) zu korrigieren. Es ergibt sich 


a oo 

e sin (S)ar+ fe °sin ar (7*) 
0 


Hier wird die Integration jedoch exakt ausführbar und führt zu dem 


Ergebnis (8%) 
— em Tolte : sin ct , coset, (wt,\\9 

1+ or} OT . 


Ist u < Tp, so folgt, wie zu erwarten, aus (8*) direkt (5*). Ist hingegen 
T < T,,so ergibt sich eine mit (9a) eng verwandte Gleichung 


R an. on 
= 


. 0o\2 2 
= ER e R (9*) 
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Damit sind auch die Gleichungen (11), (12) und (14) mit sinngemäßer 
Ersetzung von r, in (9a*) durch r, sowie A), B) und C) in der Diskussion 
durch die entsprechenden Gleichungen zu ersetzen. 


Alle anderen in der oben genannten Arbeit enthaltenen Folgerungen und 
Diskussionen bleiben von dieser Korrektur unberührt. : 


Es sei mir gestattet, an dieser Stelle besonders herzlich Herrn Dr, Van 
Vliet sowie Herrn Prof. Bittel zu danken, die mich freundlicherweise 
auf den Fehler aufmerksam machten. 


Berlin, II. Physikalisches Institut der Humboldt-Universität. 


Bei der Redaktion eingegangen am 11. August 1954 


Verantwortlich 
für die Schriftleitung: Prof. Dr. Friedrich Möglich, Berlin-Buch, Lindenberger Weg 74; 
für den Anzeigenteil: VEB Georg Thieme, Anzeigenabteilung, Leipzig C 1, Hainstr. 17-19, Aufg. C, 
Ruf 21981. Z. Z. gilt Anzeigenpreisliste Nr. 3; Verlag: Johann Ambrosius Barth, Leipzig C1, 
Salomonstraße 18B; Fernruf: 63105 und 63781. Veröffentlicht unter der Lizenznummer 285/1264 
des Amtes für Literatur und Verlagswesen der Deutschen Demokratischen Republik 
Printed in Germany Druck: Paul Dünnhaupt, Köthen (IV/5/1) I 203/54 
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ZEISS 
FUNKENERZEUGER FF 20 


nach Feussnerschem Prinzip in VDE-Aus- 
führung bei spektrochemischen Analysen 
die Gewähr für genaue Reproduzier- 


barkeit der Anregungsbedingungen 
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Anleitung zur Analyse 
der Lösungsmittel 


Von Dr. se. nat. KURT THINIUS 
Dipl.-Chemiker, Magdeburg 
86 Seiten mit 7 Abbildungen im Text. 
1953. Gr. 8°. DM 4.50 
Praktische Chemie: Diese begrüßens- 
werte Arbeit entspricht sowohl den Be- 
dürfnissen' der Gewerbehygiene wie den- 
jenigen der praktischen Chemie, zumal 
sich die Anwendungs- bzw. Geltungs- 
bereiche der Lösungsmittel-Industrien 
immer mehr ausweiten. M 


JOHANN AMBROSIUS BARTH / VERLAG / LEIPZIG 
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Die elektroakustischen Wandler q 


Von Dr. Dr. HEINRICH HECHT, Kiel 


3., erweiterte Auflage. XVI., 257 Seiten mit 52 Abbildungen i. T. 1954. Gr. 8°, 
DM 19.—, Ln. DM 20.30 
Frequenz: Dieses Buch des bekannten Verfassers gehört zum notwendigen Rüstzeug 
des Elektroakustikers. Es ist speziell dem Schallsignaltechniker »auf den Leib ge- 
schrieben«, doch kommt die allgemeine Elektroakustik deswegen nicht zu kurz; denn 
neben den meist in Resonanz arbeitenden Signalsendern und -empfängern werden auch 
der Lautsprecher, das Telephon, das Kondensatormikrophon und der Schwingquarz be- 
handelt. Deshalb gehört das Buch auf den Arbeitstisch u. =" ey: und Phy- 
sikers, der sich mit der Theorie elektroakustischer Geräte be . Gosewinkel 


Schaltschemata 
und Differentialgleichungen elektrischer und 
mechanischer Schwingungsgebilde 
Von Dr. Dr. H. HECHT, Kiel 
3., erweiterte Auflage. = 167 Seiten mit 46 Abbildungen i. T. 1954. Gr. 8°, 


M 16.—, Hin. DM 17.50 


Archiv für elektrische Übertragung: Das Buch ist nicht nur hervorragend 
geeignet, Anfänger in die Grundgedanken der technischen Schwingungslehre einzuführen, 
sondern es wird auch der erfahrene Fachmann, der bestrebt ist, sich stets immer wieder 
auch um die physikalische Fundamentierung und den systematischen Aufbau seines 
Arbeitsgebietes zu bemühen, eine Fülle von Anregungen und Wissenswertem daraus 
schöpfen. F. A. Fischer 


Die Oxydkathode 


Von Dr.-Ing. G. HERRMANN, Berlin, und Dr. phil. 8. WAGENER, London | 


I. Teil: Physikalische Grundlagen 
Von Dr. phil. S. WAGENER 
2., umgearb. Auflage. VIII, 131 S. m. 76 Abb. i. T. 1948. Gr. 8°. DM 14.40 
Vergriffen, neue Auflage in Vorbereitung 


Zeitschrift fir Naturforschung: Dieser erste Teil enthält die theoretischen 
physikalischen Grundlagen, die zum Verständnis der Oxydkathode und ihres Emissions- 
mechanismus erforderlich sind. 

In klarer Darstellung und geschickter Auswahl wird alles beriicksichtigt, was zum Ver- 
ständnis der angeschnittenen Fragen nötig ist. B.Gysae, Hechingen 


Il. Teil: Technik und Physik 
Von Dr.-Ing. G. HERRMANN, Berlin, und Dr. phil. S. WAGENER, London 


2., neubearbeitete Auflage. VIII, 284 Seiten mit 147 Abb. und 3 Tafeln. 1950. 
Gr. 8°. DM 27.—, geb. DM 29.— 


Das Elektron in Wissenschaft und Technik: Das Buch stellt eine hervor- 
g der heutigen Erkenntnisse über die Oxydkathode dar und 
kann ah oo und Physiker, der sich mit der Theorie und Anwendung von Oxyd- 
kathoden und von Verstärkerröhren beschäftigt, wärmstens empfohlen werden. Ein aus- 
führliches Literaturverzeichnis erleichtert das Einarbeiten auch in ganz spezielle Probleme 
der Kathode. F. Malsch 


JOHANN AMBROSIUS BARTH / VERLAG / LEIPZIG 
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